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ABSTRACT 


The  theory  of  virtual  cathode  oscillation  and  a  detailed  experimental 
design  of  a  millimeter  wave  virtual  cathode  oscillator  (VIRCATOR)  are  pre¬ 
sented.  It  is  found  that  second  order  analytical  approximations  give  an 
accurate  value  for  the  space-charge  limiting  current  of  an  electron  beam  in 
cylindrical  geometry.  The  same  Green's  function  approach  could  be  used  on 
other  geometries  as  well.  Exceeding  the  space-charge  limiting  current  is 
the  necessary  condition  for  virtual  electrode  formation.  However,  it  does 
not  explain  the  nature  of  the  space-charge  limit  instability.  The  oscilla¬ 
tory  stable  state  of  the  virtual  cathode  is  found  to  originate  at  the 
bifurcation  point  of  several  stable  and  unstable  electron  beam  states.  An 
analysis  of  the  nonlinear  growth  rate  of  the  instability  is  derived. 

The  virtual  cathode  is  in  effect  a  relaxation  oscillator.  The  advant¬ 
age  of  the  vircator  over  a  reflex  klystron  is  that  the  virtual  cathode 
moves  opposite  to  the  electron  beam  propagation  direction  during  the  bunch¬ 
ing  phase  resulting  in  enhanced  electron  densities  unobtainable  otherwise. 
The  periodic  time-dependent  nature  of  the  bunching  results  in  efficient 
microwave  generation.  Experimental  efficiencies  of  12%  and  theoretical 
efficiencies  of  20%  have  been  obtained.  The  microwave  frequency  scales  as 
the  relativistic  beam  plasma  frequency.  This  can  be  adjusted  via  an 
externally  applied  axial  magnetic  field.  The  result  is  an  oscillator  which 
is  frequency  agile  over  an  order  of  magnitude  without  externally  changing 
the  physical  conf i gurati on.  Output  powers  of  from  kilowatts  to  gigawatts 
are  obtainable  based  on  the  power  source.  Experimental  results  are 
reviewed  and  explained.  The  rationale  for  using  a  foilless  diode  as  the 
electron  beam  source  is  given. 

A  detailed  mechanical  and  electrical  design  is  described.  It  entails 
placing  the  vircator  in  a  high  vacuum  vessel  with  flexibility  for  further 
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development.  The  electrical  power  is  pulse  formed  using  a  lumped  artific¬ 
ial  Blumlein  line.  The  axial  magnetic  field  is  critically  damped  for  com¬ 
ponent  safety  and  circuit  design  simplicity.  The  beam  source  is  a  foilless 
diode  with  an  annular  limiter  for  collimation  and  prevention  of  electrons 
reflexing  into  the  diode.  A  W-band  (75-110  GHz)  microwave  grating  spec¬ 
trometer  has  been  designed  and  constructed  for  diagnosis  of  the  electromag¬ 
netic  wave  signal.  The  vircator  design  goal  is  1  Megawatt  at  100  GHz. 
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I.  INTRODUCTION 


Since  the  discovery  of  the  Chi ld-Langmui r  relation1*2  it  has  been 
known  that  exceeding  the  limiting  current  of  a  diode  leads  to  the  develop¬ 
ment  of  a  virtual  cathode.  Subsequently,  numerous  papers  were  written  on 
experiments  and  theory  relating  to  space  charge  limited  flows.  Reference  3 
provides  an  excellent  background  and  bibliography.  More  recently,  the 
exact  steady  state  solutions  for  electron  beams  in  one-dimensional  relativ¬ 
istic  diodes4  and  bounded  drift  spaces5  were  derived.  It  is  easily  seen 
that  for  sufficiently  large  currents  there  exist  two  steady  states  for  an 
electron  beam,  only  one  of  which  is  stable.6  At  the  space-charge  limit 
(SCL)  these  two  states  coalesce  and  above  the  SCL  they  disappear.7  As  cur¬ 
rent  is  increased  past  the  SCL,  the  beam  develops  a  jump  instability  and 
relaxes  into  an  oscillating  state. 

In  the  early  1960's,  computer  models  were  developed  which  quantita¬ 
tively  depicted  the  nonlinear  oscillatory  nature  of  the  virtual 
cathode.8"10  These  were  one-dimensional,  non-relati vistic,  electrostatic, 
multiple  sheet  models.  References  8  and  9  qualitatively  pointed  out  many 
interesting  dependencies  of  the  oscillation  frequency  and  potential  minimum 
position  on  injected  curent,  thermal  spread  and  circuit  resistance.  Refer¬ 
ence  10  presents  computer  experiments  with  one  and  two  species. 

The  phenomenon  of  virtual  cathode  formation  in  intense  relativistic 
electron  beams  figures  prominently  in  a  number  of  high  interest  research 
areas.  Devices  used  to  produce  high  current  ion  beams  for  inertial  con¬ 
finement  fusion--pinch  reflex  diode11*12  and  reflex  triodes13* 14-- 
depend  on  the  virtual  cathode  to  inhibit  electron  transport  and  use  its 
potential  well  to  accelerate  ions.  The  recent  concept  of  the  spherical 
el ectron-to-i on  converter15  requires  a  virtual  cathode. 
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The  virtual  cathode  plays  a  dominant  role  in  areas  other  than  product¬ 
ion  of  light  ion  beams  for  fusion.  It  is  attributed  with  the  main  role  in 
collective  ion  acceleration  in  neutral  gas.16*17  Control  of  virtual 
cathode  motion  is  the  mechanism  for  acceleration  in  the  Ionization  Front 
Accelerator.18*19  It  is  also  the  acceleration  method  in  two  concepts  for 
collective  effect  accelerators.20*21  A  final  application  is  the  use  of 
virtual  cathode  oscillations  to  produce  high  power,  short  wavelength  micro- 
waves.22"25  Experiments  using  reflex  triodes  have  already  been  produced 
1.4  GW  of  power  with  12%  beam-to-RF  efficiency.28  It  is  this  latter  appli¬ 
cation  which  will  be  considered  in  this  report. 

The  report  is  organized  as  follows.  Section  II  provides  an  in-depth 
review  of  the  analytical  approach  to  determining  the  space-charge  limiting 
current.  It  is  shown  by  comparison  with  particl e-i n-cel 1  simulations  that 
this  value  can  be  predicted  quite  accurately.  Section  III  analyzes  the 
nature  and  nonlinear  growth  rate  of  the  space-charge  limit  instability 
which  leads  to  the  formation  of  a  virtual  cathode.  In  Section  IV  the  basic 
considerations  in  the  application  of  the  virtual  cathode  for  generating 
coherent  microwaves  is  presented.  Finally,  Section  V  gives  a  detailed 
design,  both  mechanical  and  electrical,  of  the  millimeter  wave  vircator 
being  built  for  this  project. 


2 


II.  THE  SPACE-CHARGE  LIMITING  CURRENT 


A.  Therrv 

The  concept  of  limiting  current  can  be  visualized  in  the  following 
jver-simpl i fied  physical  picture.  As  particles  from  the  electron  beam  (see 
Figure  1)  enter  the  drift  tube,  the  charge  which  they  carry  creates  a 
potential  barrier  against  further  transmission  of  current.  If  the  injected 
current  is  greater  than  the  limiting  current,  the  potential  energy,  e<t, 
will  exceed  the  kinetic  energy  of  the  bean  (y0  -  I)mc2,  where 
e  =  electronic  charge 
<t  =  electrostatic  potential  <_  0 
y  =  relativistic  factor  (1  -  e2)-1/2 
Yo  =  value  of  y  at  injection 

6  =  v/c 

v  =  electron  velocity 
m  =  electron  rest  mass 

At  that  position  the  barrier  is  large  enough  to  stop  the  current  and  cause 
reflection  of  beam  electrons.  Thus,  the  deep  potential  well  is  called  a 
virtual  cathode.  Note,  however,  that  since  the  amplitude  of  the  well 
fluctuates  in  time  and  space,  the  beam  is  never  completely  cut  off  by  the 
virtual  cathode. 

The  general  space-charge  limiting  current  problem  does  not  lend  itself 
to  analytical  solution.  Instead  theoretical  work  on  the  subject  has 
attempted  to  solve  for  limiting  currents,  which  develop  in  REBs,  in  the 
strictly  electrostatic  case.  The  solutions  derived  are  only  valid  for 
beams  in  steady  state  equilibrium,  i.e.,  I0,  <  I|_.  Also,  the  deriva¬ 
tions  depend  on  the  following  assumptions: 
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1.  The  electon  beam  is  a  cold  plasma,  Te  =  0.  It  implies  the  elec- 
trons  are  col'.isionless,  ana  it  allows  one  to  neglect  the  vPe  and  viscos¬ 
ity  terms  in  the  momentum  transfer  equation.  Here,  Te  is  the  electron 
temperature  and  Pe  the  electron  pressure. 

2.  The  electron  beam  is  axially  homogeneous  and  infinite.  In  real 
terms  the  infinite  beam  assumption  implies  L  >  R,  where  L  is  the  drift  tube 
’ength,  and  R  is  tne  drift  tube  radius.  This  allows  the  beam  to  reach  a 
steady  state  far  enough  downstream  from  the  injection  plane  to  make  the 
equilibrium  beam  density,  n£,  velocity,  v°,  and  electric  field,  E° 
independent  of 

3.  The  equilibrium  radial  density  and  velocity  profiles  are 
azimuthally  symmetric  about  the  magnetic  axis.  In  conjunction  with 
paragraph  2  this  means  that  n°,  v°  and  Eu  are  functions  only  of  the  radial 
distance,  r. 

4.  The  electron  particle  density  and  kinetic  energy  are  independent 
of  r  at  the  anode  plane.  Therefore,  the  beam  is  uniform  in  space  and 
monoenergetic  at  injection  into  the  drift  tube--anode  foil  scattering  is 
neglected. 


5.  The  externally  imposed  longitudinal  magnetic  field,  B0,  is 

effectively  infinite.  If  this  assumption  is  not  made,  the  current-induced 

equilibrium  sel f -magnetic  fields  would  play  a  major  role  in  determining  the 

equilibrium  radial  profile  of  the  beam.  Since  B0  is  taken  as  infinite, 

the  electrons  are  guided  and  contained  by  it,  and  the  self-fields  can  be 

treated  as  perturbations.  Another  way  of  expressing  this  assumption  is  to 

♦ 

state  that  the  Larmor  radius  of  an  electron  in  the  B0  field  must  be  much 
less  than  the  beam  radius,  rj_  <  r&. 
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Since  assumptions  1  through  5  reduce  the  problem  to  the  electrostatic 
case,  solving  Poisson's  equation  will  give  the  desired  result.  Using  this 
approach  several  researchers  have  derived  expressions  for  the  space-charge 
limited  current  of  electron  beams  in  vacuum  and  in  charge  neutral  plasmas. 
Excellent  reviews  have  been  written.27*28  The  most  widely  cited  results 
are  those  of  Bogdankevich  and  Rukhadze.27  Using  the  conservation  of  energy 
and  momentum  they  are  able  to  write  Poisson's  equation  for  a  cyl i ndri cal ly 
symmetric  electron  beam  as 

0  <  r  <  r 
—  o 

r  <  r  <  R 
o  — 

with  appropriate  boundary  conditions.  Here,  j  is  the  current  density.  The 
radial  cross  section  of  the  beam  is  given  in  Figure  2a,  where  R  is  the 
drift  tube  radius,  and  r0  is  the  radius  of  the  axially  centered  beam. 

Equation  (1)  is  a  nonlinear  differential  equation.  The  approach  taken 
in  Reference  27  is  to  obtain  solutions  by  making  certain  analytical  approx¬ 
imations.  These  approximations  are  related  to  beam  energy  at  injection 
(the  nonrelati vistic,  y  =  1,  or  ultrarelativistic,  y  »  1,  limits)  and  beam 
geometry  (the  pencil  beam.  In  (R/r0)  >  1,  or  fat  beam,  £n  (R/r0)  C  1 
limits.  The  result  is  a  series  of  first  order  analytical  approximations 
which  are  restricted  in  validity  tc  narrow  regimes.  Indeed,  the  most  often 
used  relation  is  an  interpolation  formula  which  utilizes  the  uniform  dens¬ 
ity  approximation.  It  neglects  the  radial  variations  in  $,  y,  and  number 
density,  n.  This  equation  consistently  underestimates  the  true  value  of 
limiting  current. 
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'  GROUNDED  DR 1  FT  TUBE 


Figure  2.  Radial  cross  section  of  two  electron  beams. 

(a)  Solid  beam  where  r0  Is  the  beam  radius 
and  R  Is  the  drift  tube  radius,  (b)  Annular 
beam  where  a  is  the  Inner  radius,  b  the  outer 
radius,  R  the  drift  tube  radius  and 
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The  work  of  Genoni  and  Proctor29  produces  a  more  accurate,  less 
restricted  relation  for  space-charge  limiting  current.  In  particular,  it 
is  a  second  order  analytical  expression  which  takes  into  account  the  radial 
dependence  of  y,  $,  and  n.  Also,  the  beam  is  neither  confined  to  the 
ultrarelativistic  limit,  nor  is  the  beam  geometry  valid  in  only  certain 
regimes.  The  general  electron  beam  radial  configuration  is  given  in  Figure 
2b.  As  seen  from  the  figure,  R  is  the  drift  tube  radius,  a  and  b  are  the 

inner  and  outer  beam  radii,  respectively,  and  e  =  . 


Their  derivation  also  begins  with  conservation  of  energy  and  momentum, 
so  that  Poisson's  equation  for  a  cyl i ndrical ly  symmetric  beam  takes  the 
form 


1  d 

r  "37 


r 


Mrl 

dr 


!4v[(b/R)2  -  (a/R)2]  ,  . 

—  e{r)  -- 

0  elsewhere 


(2) 


subject  to  the  boundary  conditions 


y(R)  =  yq  y(0)  _>  1  and  finite  (3) 

where  v  =  Ie/mc3  is  dimensionless  current  independent  of  r  and  z.  This  is 
a  nonlinear  differential  equation,  as  expected,  showing  the  radial  depend¬ 
ence  of  y  By  choosing  an  appropriate  Green's  function  one  can  express 
equations  (2)  and  (3)  as  an  integral  equation  and  obtain  a  numerical  solu¬ 
tion.  More  importantly,  by  assuming  a  uniform  density  radial  profile  of  y 
inside  the  beam  as  a  zeroth  order  guess,  one  can  iterate  on  the  integral 
equation  to  obtain  a  second  order  analytical  approximation  to  the  current. 
Thus,  not  only  can  the  limiting  current  be  found  by  determining  the  value 
of  y  which  maximizes  the  current,  but  also  the  current  in  the  drift  tube 
for  any  value  of  y  and  r  can  be  evaluated.  Using  this  approach  an  equiva¬ 
lent  integral  equation,  valid  for  a  _<  r  _<  b,  is 
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r  (r )  =  y0  - 


u  -6 


jln  T  l  WJ 


X 

elxl 


In  —  dx 
r 


A  first  approximation  y(r)  *  yl(r)  is  obtained  by  putting  e  =  K  (a  con¬ 
stant)  in  the  integrals  of  equation  (4).  If  K  is  fixed  by  requiring 
y*(a)  =  (1-K2)-1/2,  one  obtains  the  following  approximation  to  the  lim¬ 
iting  current: 


where 


In  (1  -  e' 


Equation  (5),  which  appears  in  a  paper  by  Miller  and  Straw30,  is  the  logi¬ 
cal  general ization  to  annular  beams  of  the  widely  used  Bogdankevich- 
Rukhadze27  interpolation  formula  for  solid  beams,  to  which  it  reduces  in 
the  limit  e  -*■  1.  A  second  approximation,  yU(r),  is  obtained  by  repeat¬ 
ing  this  procedure  using  B  =  [1  -  (y*)-2]1/2  in  the  integrals  of  equa¬ 
tion  (4).  In  the  special  case  of  a  solid  beam  (a  =  0),  the  integrals  can 
be  evaluated  in  closed  form,  but  for  the  general  case  of  an  annular  beam 
they  must  be  done  numerically.  Genoni  and  Proctor29  have  obtained  useful 
analytical  approximations  to  yll(r)  annular  beams  by  doing  the  second 
iteration  with  b  =  (1  -  f2)1/2  where  r(r)  is  a  quadratic  fit  to  y^r) 
which  yields  the  exact  y^(r)  in  the  solid  beam  limit.  A  particular 
choice  of  K  then  yields  a  definite  relation  between  v  and  ya  =  y(a).  The 
simplest  of  these  that  gives  good  agreement  with  numerical  results  is  evi¬ 
dently 


u?  -  v1'*  -  k2  ■ 


(1  -  f  +  2  In  R/b)‘ 


9 


in  which 


r 


b 


Ya  +  •  V 


T~~f  TT  h  R/b 


(8) 


and 


9(e) 


In 


1  -  e/2 
"1  -  e 


(9) 


In  tne  limit  t  *  I,  both  g  and  f  tend  to  zero,  and  equation  (7)  reduces  to 
a  formula  (equation  14)  obtained  by  Thompson  and  Sloan  (Ref.  8)  in  their 
analysis  of  solid  beams.  It  should  be  noted  that  different  quadratic  fits, 
while  producing  more  complicated  analytical  expressions,  produce  more 
accurate  results.29 


B.  Simulation 

The  space-charge  flow  simulations  were  carried  out  using  a  two-dimen¬ 
sional  particle-in-cell  plasma  simulation  code,  CCUBE,  which  is  fully  rela¬ 
tivistic  and  electromagnetic.32*33  CCUBE  (Version  One)  solves  self- 
consistently  for  the  time-dependent  trajectories  of  tens  of  thousands  of 
plasma  particles  over  thousands  of  plasma  periods.  The  critical  portion  of 
the  code  consists  of  a  leap  frog  procedure.  At  each  time  step,  the  charge 
and  current  associated  with  the  particles  are  interpolated  bilinearly  onto 
a  spatial  mesh.  The  quantities  then  serve  as  source  terms  in  Maxwell's 
equations  for  determining  the  electromagnetic  scaler  and  vector  potentials 
on  the  mesh.  Finally,  the  potentials  are  interpolated  back  to  the  particle 
positions.  The  basic  coordinate  system  is  cylindrical  with  azimuthal  sym¬ 
metry  assumed.  A  uniform  axially  directed  magnetic  field  can  be  imposed  to 
insure  beam  equilibrium. 

All  variables  are  expressed  in  dimensionless  terms.  In  particular, 
the  depth  of  the  potential  well  formed  by  the  electronic  space-charge  is 
given  by  4,  where 
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(10) 


♦  (f)  =  =  t(0  -  T0  <  0 

me 

The  five  components  in  configuration-momentum  space  are  XI  =  z  /c,  X2  = 
r  Up/c,  VI  =  y  $2 ,  V2  =  y  0p,  V3  =  y  Be  where  Bz»Br»  and  eQ  are  the  compon¬ 
ents  of  0  in  the  z,  r,  and  e  directions,  respectively,  w®  =  (^np11^2^)1^ 
is  the  electron  plasma  frequency  arising  from  the  beam  density  at 
injection,  n^. 

In  these  simulations  the  code  is  made  to  model  the  injection  of  cold, 
monoenergetic,  relativistic  electrons  with  uniform  density  into  a  long, 
evacuated,  equipotential ,  cylindrical  drift  tube  along  an  effectively  infi¬ 
nite  magnetic  guide  field  in  order  to  comply  with  the  assumptions  of  the 
previous  section.  This  is  accomplished  in  the  following  manner. 

1.  A  cold  bean  is  obtained  by  setting  the  electron  thermal  velocity 
components  equal  to  zero. 

2.  Miller  and  Straw34  provide  a  criterion  for  neglecting  the  presence 
of  end  plates  and  effectively  having  L  »  R. 

L/R  >  2.58  (b/R)° -133  (11) 

For  the  annular  beam  simulations  this  yields  L/R  2.50,  whereas  the  worst 
case  simulation  is  L/R  =  9.84.  In  addition,  two  runs  were  made  with  the 
downstream  axial  boundary  grounded.  The  presence  of  end  plates  did  not 
affect  the  results. 

3.  Azimuthal  symmetry  is  assured  by  the  nature  of  the  code. 

4.  The  code  injection  can  be  made  to  neglect  anode  foil  scattering, 
so  that  the  particles  are  monoenergetic  and  spatially  uniform  at  inject¬ 
ion. 

5.  Thode,  et  al.35,  give  a  numerical  relation  to  insure  that  the 

imposed  magnetic  field  is  effectively  infinite.  It  is  w  /w°  >  5  where 

c  p 

oac  is  the  beam  cyclotron  frequency  defined  as  eB0/mc.  This  is  met  by 
running  the  code's  dimensionless,  external,  magnetic  field  as 
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The  number  of  cells  in  the  XI  and  X2  direction  were  190  and  37, 
respectively.  The  XI  length  was  held  constant  at  50.0  units  giving  aXI  = 
0.263  per  cell.  The  X2  radius  varied  from  2.89  to  5.0  depending  on  the 
value  of  injected  current.  The  boundary  conditions  for  the  electrons  are 
reflection  at  the  axis  and  absorption  on  both  the  radial  and  axial  sur¬ 
faces.  Potentials  are  set  to  zero  at  the  upstream  axial  and  outer  radial 
boundaries  in  order  to  represent  metal  surfaces.  At  the  downstream  axial 
boundary,  however,  the  normal  derivatives  of  the  potentials  are  set  to  zero 
in  order  to  mock  up  an  infinitely  long  cylinder. 

In  the  next  section,  the  value  of  ya  from  the  simulation  calculation 
is  compared,  for  a  given  value  of  v,  with  that  given  by  equation  (7),  as 
well  as  with  a  numerical  solution  obtained  by  iterating  on  equation  (4) 
with  r  =  a.29 

C.  Results 

Several  CCUBE  runs  were  made  using  y0  =  5  and  varying  values  of  v  in 

order  to  determine  the  maximum  value  of  |$|.  This  occurred  consistently  at 

the  radial  value  r  =  a.  Therefore,  |<t>|  maximizes  for  $(a)  =  $a.  The 
characteristics  of  the  simulation  runs  are  given  in  Table  1.  The  last  col¬ 
umn  in  Table  1  refers  to  whether  the  downstream  axial  boundary  was  grounded 
or  open.  Such  a  condition  enabled  us  to  check  whether  or  not  end  plate 
effects  were  entering  the  problem.  As  seen  in  Table  1  no  effects  were 
noted.  The  percent  error  was  derived  on  the  basis  of  the  code's  management 
of  conservation  of  energy,  which  is  denoted  by  a  net  energy  diagnostic. 
Since  net  energy  continuously  increases  due  to  roundoff  and  truncation 
errors,  the  percent  error  will  increase  as  a  function  of  computer  run 
time.  The  percent  errors  given  in  Table  1  were  taken  after  the  run  reached 

a  steady  state,  as  denoted  by  total  energy  and  |$|  reaching  a  constant,  or 

a  virtual  cathode  was  formed.  The  latter  case  was  true  in  runs  3  and  8. 


This  can  be  observed  in  the  table,  since  |<f>aj  >  (yQ  -  1).  Thus,  the 
space-charge  limiting  current  lies  between  v  =  1.175  and  1.293  for  the 
solid  beam  case  and  v  =  2.846  and  3.034  for  the  annular  beam  case. 

TABLE  1 

CHARACTERISTICS  OF  SIMULATION  RUNS 


Downstream 


Run 

V 

hi 

Error 

(%) 

a/R 

b/R 

L(c/u>°) 

R(c/uj°)  Axial 
p  Boundary 

1 

1.175 

3.07 

1.5 

0 

0.5 

50.0 

4.38 

0 

2 

1.175 

3.06 

1.2 

0 

0.5 

50.0 

4.38 

G 

3 

1.293 

4.20 

2.5 

0 

0.5 

50.0 

4.60 

0 

4 

0.982 

1.01 

0.74 

0.4 

0.8 

50.0 

2.89 

0 

5 

1.957 

2.05 

1.1 

0.4 

0.8 

50.0 

4.08 

0 

6 

2.846 

3.18 

1.2 

0.4 

0.8 

50.0 

4.92 

0 

7 

2.846 

3.18 

1.1 

0.4 

0.8 

50.0 

4.92 

G 

n 

o 

3.034 

4.25 

3.1 

0.4 

0.8 

50.0 

5.08 

0 

A 

comparison  is  made 

between 

the  numerical, 

analytical  and  sinualtion 

models  in  Table  2.  All  three  methods  produce  consistent  values  of  ya  for 
given  values  of  sublimiting  injected  current.  Table  2  lists  the  values 
given  by  the  three  methods,  where 


is  used  tc  define  ya  for  the  simualtions.  Note  that  if  a  virtual  cathode 
is  formed,  j$a|  >  {y0  -  1),  ya  as  defined  by  equation  (13)  does  not 
have  a  physical  meaning.  Rather,  it  is  simply  a  mathematical  relation 
between  y0  and  j<j>a|  • 
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TABLE  2 

COMPARISON  OF  MODELS 


V 

Ya 

Numerical 

Ya 

Analytical 

Ya 

Simulation 

Beam  Type 

1.175 

1.893 

1.93 

1.93 

Solid 

1.293 

Vi rtual 
Cathode 

Vi  rtual 
Cathode 

Vi rtual 

Cathode 

Solid 

0.9  82 

4.004 

4.00 

3.99 

Annular 

1.957 

2.976 

2.98 

2.95 

Annular 

2.846 

1.867 

1.92 

1.82 

Annular 

3.034 

Vi rtual 
Cathode 

Vi  rtual 
Cathode 

Vi  rtual 

Cathode 

Annular 

A  more  important  comparison  can  be  made  between  the  values  for  space- 
charge  limiting  current,  v|_,  derived  from  the  numerical,  analytical,  and 
simulation  models  and  that  produced  by  the  widely  used  Bogdankevich- 
Rukhadze  interpolation  formula,  equation  (5),  This  is  done  in  Table  3. 
Equation  (5)  predicts  values  for  vl  which  are  significantly  lower  than 
those  of  equation  (4)  and  (7),  However,  the  analytical  results  differ  from 
the  figures  given  by  the  integral  equation  by  roughly  3  percent.  v|_  is 
obtained  from  equation  (4)  by  determining  numerically  the  largest  value  of 
v  for  which  a  solution  exists.  The  limiting  current  for  equation  (7)  is 
derived  from  the  value  of  ya  which  maximizes  v. 

TABLE  3 

.  LIMITING  CURRENT  COMPARISON 


VL 

VL 

VL 

VL 

Beam  Type 

Numerical 

Analytical 

Simulation 

B-R 

Solid 

1.27 

1.31 

1.175  -  1.293 

1.118 

Annular 

3.03 

3.12 

2.846  -  3.034 

2.712 

The  discrepancy  in  the  limiting  value  of  v  between  the  Bogdankevich- 
Rukhadze  and  Genoni -Proctor  formulations  is  due  to  the  error  involved  in 
assuming  the  uniform  density  approximation.  Figure  3  is  a  plot  of  the 
relation  between  v  and  ya  in  the  uniform  density  approximation  and  in  a 
second  order  analytical  expression  using  a  better  quadratic  fit  than  equa¬ 
tion  (7).  The  simulation  points  are  added.  It  is  clear  that  the  uniform 
density  approximation  underestimates  the  value  of  v  for  all  values  of 
ya.  Finally,  Figure  4  plots  y(r)  -  1  versus  r.  It  shows  the  dependence 
of  beam  kinetic  energy  on  r  for  the  annular  beam  modeled.  The  successively 
more  energetic  lines  of  dots  correspond  to  the  more  recent  injections  of 
charge  through  anode.  As  expected,  the  minimum  kinetic  energy  occurs  at 
the  inner  beam  edge,  and  the  maximum  occurs  on  the  outer  edge.  The  unifrom 
density  approximation  neglects  this  dependence. 
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Figure  3.  Plot  of  v  versus  y»  for  y0  *  5,  a/R  »  0.4,  b/R  *  0.8.  The  dashed  line 
represents  the  uniform  density  approximation.  The  solid  line  results 
from  an  approximate  second  order  analytical  relation.  The  points 
correspond  to  simulation  results.  (N.B.  ya  for  the  point  labeled  v.c. 
(virtual  cathode) has  no  physical  meaning.  It  is  defined  in  terms  of 
equation  (13)). 


w*p  t=260 . 0800 

Figure  4.  A  CCUBE  diagnostic  of  particle  kinetic  energy 
versus  r  for  the  annular  beam  showing  the 
radial  dependence  of  y. 
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III.  THE  SPACE-CHARGE  LIMIT  INSTABILITY 


We  have  visualized  the  space-charge  limiting  current  as  a  point  at 
which  the  electrostatic  potential  exceeds  the  kinetic  energy  of  the  beam. 
However,  it  is  easily  noted  from  Figure  3  that  at  the  SCL  the  beam  kinetic 
energy  is  approximately  y01/3mc2.  Thus,  although  the  onset  of  virtual 
cathode  formation  can  be  exactly  predicted  as  the  point  at  which  the  beam 
exceeds  a  critical  current,  the  dynamics  of  its  formation  and  subsequent 
oscillations  are  not  well  understood. 

A  small  signal  perturbations  analysis  can  be  performed  just  below  the 
SCL  which  seems  to  predict  instability  above  the  limiting  current 
value,6*9  however,  it  is  not  rigorous  and  in  fact  produces  misleading 
results.  Instead  one  can  use  multiple  scaling  perturbation  techniques  to 
study  the  time  dependent  behavior  of  a  beam  when  the  SCL  is  exceeded.  We 
derive  estimates  for  the  nonlinear  growth  rate  of  the  ensuing  instability 
and  show  that  even  below  the  SCL  the  beam  is  unstable  to  sufficiently  large 
perturbations.  The  method  can  be  applied  to  a  wide  class  of  problems,  but 
here  we  treat  the  short-circuited  one-dimensional  electrostatic  diode 
depicted  in  Fig.  5  as  the  simplest  model  containing  the  appropriate  phys¬ 
ics.  We  show  that,  at  least  in  one  dimension,  an  arbitrarily  heavy  ion 
background  does  not  alter  the  qualitative  behavior  of  the  beam,  and  present 
numerical  results  that  exhibit  virtual  cathode  oscillations  for  a  neutral 
beam.  Before  starting  this  analysis,  it  is  worthwhile  to  describe  the 
results  of  numerical  simulations  to  present  a  picture  of  time-dependent 
virtual  cathode  dynamics. 
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A.  Physical  Description  of  the  Virtudl  Cathode 


Simulations  were  carried  out  in  conjunction  with  the  theory  presented 
in  the  next  section  using  a  two-dimensional,  relati vistic,  electrostatic, 
particle-i n-cell  code.  The  code  can  solve  sel f-consi stently  for  the  time 
dependent  trajectories  of  tens  of  thousands  of  plasma  particles  over  thou¬ 
sands  of  plasma  periods.  All  variables  are  expressed  in  dimensionless 
terms.  Therefore,  length  is  in  units  of  c /ujp ;  time  is  measured  in  units 
of  u>p~l,  and  particle  velocity  is  given  by  v-j  =  Bjy  (i  =  1,2,3), 
where  up  is  the  initial  electron  plasma  frequency. 

In  these  simulations  a  monoenergetic  51  keV  electron  beam  is  injected 
into  a  Cartesian  geometry.  The  left  and  right  boundaries  are  grounded  rep¬ 
resenting  a  planar  short-circuited  diode.  Periodic  boundary  conditions  in 
the  transverse  direction  make  configuration  space  effectively  one-dimen¬ 
sional.  In  general,  the  simulation  had  62  cells  in  the  longitudinal 
direction  modelling  a  length  of  1.0  c/wp.  The  time  step  was  0.0125 
Twenty  particles  were  injected  per  cell. 

A  detailed  discussion  of  the  physical  dynamics  of  the  virtual  cathode 
based  on  these  numerical  results  is  appropriate  here.  The  usual  graph  of 
potential  minimum,  <f>m,  in  the  diode  versus  electron  beam  current,  o,  is 
shown  in  Fig.  6.  The  arameter  a  will  be  discussed  later.  When  a  is 
increased  above  the  space-charge  limit,  jumps  from  the  stable  normal-C 
branch  to  the  oscillatory  stable  branch.  The  amplitude  and  position  of 
while  on  the  oscillatory  branch  describe  a  limit  cycle,  as  expected 
for  a  relaxation  oscillation  which  this  represents.  Typical  limit  cycles 
are  depicted  in  Fig.  7.  As  a  is  increased  further,  <f>m,  the  oscillation 
frequency,  and  virtual  cathode  position  within  the  diode  asymptotically 
approach  limiting  values.  If  a  is  decreased,  the  oscillation  amplitude, 

A4>m ,  decreases  and  the  position  of  moves  toward  the  diode  center. 
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NORMAL-C 


(stable).  The  motion  around  the  hysteresis  loop  is  denoted  by 


51 

is 


The  electron  flow  reverts  to  the  equilibrium  steady  state  when  the  pertur¬ 
bation  due  to  the  rate  of  change  of  diode  current  below  the  space-charge 
limit  is  sufficiently  large.  This  normally  occurs  before  the  bifurcation 
point  is  reached.  The  entire  process  forms  a  hysteresis  loop,  which  is 
depicted  in  Fig.  6. 

The  virtual  cathode  originates  at  the  bifurcation  point.  This  is  the 
intersection  of  the  oscillatory  state  with  the  C-overlap7  and  partially 
reflected  solution  branches.5  The  bifurcation  point  cannot  be  reached  in 
the  short-circuited  diode.  Of  the  three  branches  emanating  from  it  two 
(the  steady  ones)  are  physically  unstable  while  the  oscillatory  branch  is 
numerically  unstable  at  this  point.  This  results  because  the  limit  cycle 
at  the  bifurcation  point  is  infinitesimally  small,  so  that  simulation  codes 
lose  resolution  before  it  can  be  reached.  Loss  of  resolution  creates  a 
small  amplitude,  high  frequency  oscillation  observed  in  this  study  and  pre¬ 
viously.9  This  result  is  numerical,  not  physical. 

This  problem  can  be  overcome,  if  we  eliminate  the  hysteresis  loop. 

Then  the  C-overlap  branch  disappears  and  we  can  get  to  the  bifurcation 
point  along  the  stable  normal -C  branch.  This  can  be  accomplished  in  sev¬ 
eral  ways.  The  most  appropriate  in  this  study  is  to  have  a  retarding 
potential  difference  across  the  diode  equal  to  the  injected  electron 
kinetic-energy.  Then  the  C-overlap  solution  vanishes  and  the  bifurcation 
takes  place  at  the  space-charge  limit,  a^.  The  oscillation  can  be 
described  as  a  small  perturbation  on  the  beam  rather  than  the  radical 
change  that  results  in  the  short-circuited  diode  when  aS£L  is 
exceeded.  This  analysis  indicates  that  the  onset  of  virtual  cathode  forma¬ 
tion  occurs  when  the  electron  velocity  in  the  steady  state  vanishes  at  some 
point  inside  the  diode.  For  the  short  circuited  diode  this  occurs  at  the 
diode  center;  for  the  biased  diode  it  occurs  at  x  <  £.  Because  the  oscil¬ 
lation  is  a  small  perturbation  on  the  steady  state  fields  in  the  biased 
diode,  it  is  readily  seen  that  the  virtual  cathode  oscillation  period  at 


onset  is  the  electron  transit  time  from  the  injection  plane  to  the  position 
where  velocity  vanishes. 


Finally,  consider  the  particle  dynamics  during  the  oscillation  period 


for  counterclockwise  limit  cycles 


(o  <  ajClJ  as  in  ^9*  ^c* 


At  the 


point  where  the  virtual  cathode  position  is  a  minimum  and  the  potential 


well  is  starting  to  move  to  the  right,  its  amplitude  is  too  small  to  stop 
the  electrons.  When  its  motion  is  to  the  left  it  opposes  the  electron  beam 
and  causes  particle  bunching.  Since  the  well  is  deeper,  the  stream  veloc¬ 
ity  will  vanish  at  some  location  and  then  become  negative.  Here,  the 
second  derivative  of  the  velocity  (d2u/dx2)  is  also  negative.  In  this  pro¬ 
cess  the  stream  is  continuously  deformed  to  create  a  double  valued  negative 
velocity  protrusion.  The  entire  system  is  three  valued  (Fig.  8),  as  in  a 
collisionless  shock  wave.36  Here,  the  region  of  triple  flow  is  not  limi¬ 
ted,  as  in  usual  collisionless  shocks,  by  the  presence  of  a  transverse  mag¬ 
netic  field,37  but  rather  by  the  presence  of  the  walls.  Indeed,  the 
reflected  part  detaches  from  the  main  beam  and  exits  through  the  anode  per¬ 
iodically,  thus  being  responsible  for  the  onset  of  oscillatory  behavior  in 
beam  characteristics.  As  the  potential  minimum  reaches  the  end  of  its  left 
motion  the  two  "lips"  of  the  back  reflected  stream  close.  At  this  point  no 
more  electrons  are  reflected  and  the  well  moves  to  the  right  repeating  the 


cycle. 


For  larger  values  of  a  the  limit  cycle  is  distorted  into  a  figure  8 
with  one  lobe  having  clockwise  and  the  other  counterclockwise  motion  (Fig. 
7b).  This  transition  continues  until  the  motion  is  completely  clockwise 
(Fig.  7a).  It  indicates  a  change  in  the  particle  bunching  process  and  is 
related  to  the  fraction  of  current  which  is  reflected  versus  transmitted 
from  the  injected  electron  beam.  In  Figure  7a  most  of  the  beam  is 
reflected  whereas  in  Fig.  7c  most  electrons  are  transmitted. 
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x(c/Wp) 


Figure  8.  Successive  snapshots  of  electron  bean  nomentum  space  and  corres¬ 
ponding  potential  shape  In  the  diode  for  a  =  2.0,  I  =  1.0 
c/wp.  The  time  between  frames  is  .5  wp_1.  The  initial  bean 
kinetic  energy  is  51  keV. 


e6/mc  pxy  e<t>/mc 


k  0.00 


0.00  0.25  0.60 

X  (c/w«) 


-0.05 


-0.1 


-0.150 


-0.2 


!  0.00 


X  (C/WB) 


”0.050 


-0.100 


-0.150 


-0.200  ‘ - 1 - 1 - * - - 

0.00  0.25  0.50  0.75  1.00 


Figure  8.  Concluded. 
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B.  Theory 


The  one-dimensional  motion  of  electrons  in  the  diode  is  given  by  the 
equations  of  continuity,  and  momentum  conservation  for  the  electrons  plus 
Poisson's  equation.  They  are  expressed  here  as 

pt  +  (p  v>x  =  0 

vt  +  v  vx  =  ‘  n  5 
eo  5x  =  p  +  ph 

where  0  <_  x  The  subscripts  denote  differentiation  with  respect  to 

that  variable.  Electron  and  heavy  ion  charge  density,  electron  velocity, 
electric  field,  time  and  position  are  indicated  by  p,  p^,  v,  £,  t,  and  x, 

respectively.  The  appropriate  boundary  conditions  are  v(o,t)  =  vQ, 

p(o,t)  *  pQ  and  /  5  dx  =0.  An  ion  component  is  placed  in  Poisson's 

equation  in  order  to  discuss  two-specie  space-charge  flow.  Conservation 
equations  of  mass  and  momentum  for  the  heavy  ions  are  not  included,  because 
it  is  assumed  that  their  velocity  does  not  change  appreciably  during  their 
passage  through  the  diode. 

In  order  to  simplify  the  mathematics,  we  introduce  dimensionless  equa¬ 
tions  for  our  model.  They  are 


+  (nu)x  =  0 

(15a) 

+  uux  =  -  ctE 

(15b) 

=  a(n  +  I) 

(15c) 

(14a) 

(14b) 

(14c) 
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where  n  =  p/|p0| •  u  =  v/vQ,  E  =  (eeQ/m| pQ| )l ^  C/vQ,  t  =  vQt /£,  x  =  x/i, 

*  =  <e|po|/eom>1/2  */vo  a  dimensionless  parameter  related  to  current, 

and  I  the  ratio  of  heavy  ion  to  electron  charge  densities,  I  =  Ph/ | p0 | • 

Alternatively,  a  may  be  written  as  £up/v0  where  u p  is  the  beam  plasma 

frequency.  The  boundary  conditions  for  electrons  become  u(o,  t)  =  1, 

1 

n  (o,  t)  =  -1,  and  ^  E  dx  =  0.  For  an  ion  beam  n  (o,  t)  =1,  otherwise  the 
following  derivation  is  the  same. 

1.  Unneutralized  Beams 


For  an  unneutralized  beam,  setting  I  =  0  and  solving  by  the  method  of 
characteristics38*39  we  find 

i  2  ? 

n  =  -  (t-s)  +  a  Eq(s)  (t-s)  -  1  (16) 

where  s  is  the  entry  time  for  the  particle  occupying  position  x  at  time  t 
and  E0(s)  is  the  electric  field  at  x  =  0.  The  particle  trajectories  are 
found  by  utilizing  equation  15a,  from  which  it  follows  that 

(!i)t  <17> 

This  yields 


x 


'  (t-s)  +  a  / 


EQ(s )  (t-s)  ds  +  (t-s) 


(18) 
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Integration  of  the  trajectory  equation  is  hard  for  general  time- 
dependent  situations,  because  imposing  the  proper  boundary  conditions  leads 
to  a  nonlinear  integral  equation  for  E0(t).  However,  several  special 
cases  can  be  solved  exactly.  The  problem  of  injection  into  an  empty  diode 
can  be  integrated  until  the  formation  of  a  singularity  in  n,  indicating  the 
crossing  of  trajectories.3 9  In  this  case  the  stream  velocity  becomes  3- 
valued  and  one  must  use  a  Vlasov  equation  description,40  rather  than  system 
(15)  that  is  derived  assuming  a  single  stream  of  monoenergetic  particles. 

As  described  in  the  previous  section  this  multistreaming  is  characteristic 
of  the  oscillatory  state  created  when  a  exceeds  its  SCL  value. 

Using  these  equations  we  can  derive  a  similar  representation  for  Fig. 

6  in  terms  of  E0  and  a.  For  steady  states,  E0(t)  =  E0,  a  constant, 
we  find 

u  =  -n_1  =  (t-s)*  -  cE0(t-s)  +  1  (19a) 

x  =  (t-s)3  -  f  E0(t-s)2  +  (t-s)  (ISt) 

Imposing  the  conditions  x  =  1,  u  =  1  at  t-s  =  t0,  the  particle  transit 
time,  we  note  that  t0  must  satisfy 

f?  lo  -  ‘o  +  1  =  0  (20) 

This  equation  has  two  positive  solutions  for  0  _<  a  _<  4/3,  coalescing  at 
a  =  4/3.  The  largest  one,  for  0  £  a  <  2/2/3  does  not  correspond  to  a  real 
flow.  In  Fig.  9  we  show  E0  (=  at0/2)  vs.  a.  This  representation  will 
be  used  in  the  discussion  of  nonlinear  stability. 
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2.  Neutralized  Beams 


The  Pierce  instabiity  occurs  when  electron  and  ion  space-charge  flow 
is  considered  in  finite  geometries  where  there  is  no  potential  difference 
across  the  boundaries.4 1  The  ions  can  be  stationary  or  moving  with  respect 
to  electrons.  Charge  neutrality  is  maintained  at  the  injection  plane. 

This  instabiity  may  have  ramifications  for  charged  particle  inertial  con¬ 
finement  fusion  because  of  its  effect  on  neutralized  beam  propagation  in 
the  reactor.42  In  this  section  we  show  that  the  Pierce  instability  is  a 
special  case  of  electron  and  ion  space-charge  flow.  In  general,  two 
species  flow  has  steady  and  oscillatory  states  analogous  to  one  species 
space-charge  flow. 

The  steady  state  behavior  for  the  case  of  arbitrary  I  can  be  found  in 
a  manner  similar  to  I  =  0.  Rewriting  system  (15)  in  characteristic  coor¬ 
dinates  we  arrive  at 

2 

U  +  t  )  *  C2I  (1  +  i)  -  0  (21) 

dt 

For  positive  ions  (I  >  0)  the  solution  '  (21)  after  satisfying  the  bound¬ 
ary  conditions  is 

1  +  -j-  =  (1  +  I)cosot/l(t-s)  +  Eo/Tsina/T(t-s)  (22) 

Imposing  conditions  x  =  1,  u  =  1  at  t-s  =  t0,  we  find  the  system 
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1  -  I 

~^r~ 


COSa/T  t 


(23b) 


1  = 


E 

a  "  T’  Sina/T  to 


For  1=1,  which  implies  charge  neutralization,  these  equations  redu-e 
to  the  relations  given  in  Reference  43  for  the  Pierce  instability.  How¬ 
ever,  by  varying  I  the  curves  shown  in  Figure  10  are  obtained.  These  are 
cuts  at  constant  I  through  a  three-dimensional  contiguous  surface.  The 

space  is  defined  by  the  axes  E  =  E0 1 1 /2 ,  A  =  aI3/2  and  I.  The  sur¬ 
face  is  2it  periodic  in  A  with  the  vertical  plane  at  A  =  2n  being  common  for 
all  values  of  I.  For  given  A  a  linearized  analysis  establishes  that  the 
equilibria  denoted  by  the  curves  are  stable  (unstable)  for  the  lowest 

(highest)  value  of  E.  At  1=1  exchange  of  stability  takes  place  at  odd 
multiples  of  it.  For  I  <  1  exchange  of  stability  occurs  at  the  points  where 
dE0/da  -*•  «. 


It  is  evident  from  Fig.  10  that,  for  I  <  1,  there  are  no  stable  equi¬ 
librium  solutions  in  the  neighborhood  of  A  =  w.  Therefore,  one  expects  a 
virtual  cathode  to  form  when  I  <  1  and  A  adiabatical ly  increases  to  n.  We 
have  found,  by  using  numerical  simulation,  that  in  this  case  the  beam 
settles  to  an  oscillatory  state,  similar  to  the  virtual  cathode  for  unneu¬ 
tralized  beams.40  By  slowly  increasing  I  past  the  neutral  beam  value  of  1 
in  our  simulation,  we  have  established  that  this  oscillation  persists. 
Indeed,  finding  this  oscillatory  state  for  I  >  1  by  other  means  would  have 
been  difficult,  because  the  simulation  would  tend  to  follow  the  stable 
steady  branch  that  is  present  for  all  values  of  current. 


3.  Nonlinear  Stability  Analysis 

For  I  <  1,  it  is  of  interest  to  establish  the  properties  of  the  beam 
Instability  at  the  SCI  -  generalized  for  I  *  0  to  mean  the  point  where 
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da/d£0  =  0.  We  shall  carry  out  the  analysis  for  I  =  b,  but  our  method  can 
be  applied  to  any  similar  jump  phenomenon. 


A  linear  stability  analysis6  about  the  steady  state  described  by  (20) 
results  in  the  dispersion  relation 


(2  +  8)  e"0  =  2  -  g  +  3 

a  lo 


(24) 


where  g  =  iu>t0.  We  have  written  the  expression  derived  in  Reference  6  in 
terms  of  our  dimensionless  variables.  For  a  near  the  SCL  value  we  let 


e  <<  1. 


(25) 


and  find  from  (20)  that  near  this  value,  t0  is  approximately 


3  -  3  .  cm 

7  *  e17f  0(e  1 


(26) 


where  the  -(  +  )  sign  corresponds  to  the  lower  (upper)  branch  in  Fig.  9. 

By  substituting  in  (24),  and  assuming  g  small,  we  find  that 

in)  =  g  =  +  2/T  t  +  0(e2)  (27) 

Since  the  linearized  analysis  led  to  time  factors  of  the  form  e1u)t:  in 
the  perturbations,  it  follows  that  the  lower  branch  in  Fig.  9  is  stable  and 
the  upper  unstable,  while  at  the  SCL  (e  =  0)  we  have  neutral  stability. 

Above  the  value  a  =  4/3,  linearized  theory  is  not  applicable.  Utiliz¬ 
ing  multiple  scaling  perturbation  theory41*  we  can  carry  out  a  nonlinear 


33 


stability  analysis  near  a  =  4/3.  In  system  (15)  we  set  a  =  4/3  ±  e2.  In 
this  neighborhood  perturbations  evolve  on  a  "slow"  time  scale,  depicted  by 
x  =  et. 

Eliminating  the  electric  field  by  combining  (15b)  and  (15c)  and 
utilizing  t,  system  (15)  becomes 

en  +  (nu)  =  0  (28a) 

t  '  x 

(euT  +  uux)x  =  ~  (7  1  c2)2f1  (28b) 


1 

with  conditions  u(o,t)  =  1,  n(o,i)  =  -1  and  /  Edx  =  0  rewritten  as 

o 


cju  dx  +  l  [u2(1,t)  -  u2(0,t)]  =  0 
o  T  2 

Substituting  the  asymptotic  expansions 

u  -  l  e1 u. (x,t )  +  0(eJ  +  1),  n  .  £  e*n. (x,x)  +  0(e^  +  1) 

1-0  1  i-0  1 


(29) 


(30) 


for  u  and  n  into  (28)  and  equating  coefficients  of  various  powers  of  e, 
there  results  a  hierachy  of  equations  for  the  u-j  and  n-j. 


Solving  the  0(1)  system  gives 

(u0  -£)(u0  *  l)2  «  2(2*  -  1): 


n=  -  1/u 
o  o 


(31a) 

(31b) 
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To  solve  the  0(e)  system  we  introduce  a  new  variable  q  by 


x 

so  that 


We  then  find 
nl 

U1 


T~  (§“  "  "5  q2)  +  q 

(32) 

- 

(3" ) 

±7  ul 
uo 

(34a) 

c  ala  -  Ull 

(34b) 

with  C  a  constant  of  integration  which  is,  in  general,  a  function  of  the 

slow  time  t.  To  find  C  =  C(t)  which  determines  the  slow  evolution  of  the 

2 

perturbation  u^x.t)  we  need  to  go  to  the  next  order,  0(e  ).  By 

substituting  in  the  expressions  for  uq,  nQ,  n^  and  eliminating  n^.  we 
find  that  u2(x,t)  satisfies 


(V2>xx  “2 


,  dq  ♦  c2 


0  0  u 


c  hirim 

T  l  u 


+  4  - 

0  /x  0 


(35a) 


u2 (0»t  )  =  0,  u2(1,t)  =  +  75-  CT 


(35b) 
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The  solution  to  this  inhomogeneous  two-point  boundary  value  problem  exists 
provided  a  certain  orthogonality  condition  between  the  right  hand  side  and 
the  solution  of  the  adjoint  problem  that  takes  account  of  the  boundary  con¬ 
ditions  is  satisfied  (Fredholm  alternative  theorem).45  This  leads  to  the 
desired  equation  determining  C(t): 

aCT  +  bC2  t  c  =  0  (36) 


where  a,  b,  and  c  are  found  to  be 


=  .  j/z  q(q-3/2) 

0  U 


^q2(q  -  9/4)  +| 


dq  =  1.6850 


b  *  |  },Z  (13(^~3./2I3  dq  =  -  3.7968 

0 


3/2 

c  *  4  /  q(q-3/2)dq  =  -1.5 

J  o 


(37a) 


(37b) 


(37c) 


In  (36)  the  (+)  or  (-)  signs  indicate  that  we  are  above  or  below  the  SCL, 
respecti vely. 


Above  the  SCL,  we  find 


) 


(38) 


and  below 


C(0  -  F 


tanh 


/cF  (^)  if  | C ( o ) | 


<  /c/b 


(39) 
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|C(o)|  >  /c/b 


(40) 


C(t) 


coth 


if 


where  t0  is  a  constant  of  integration.  In  general,  small  initial  pertur¬ 
bations  will  lead  to  the  solution 


u(x,t) 


uo+I  Ci(T)et°it  ui(x)+0(e2) 


(41) 


where  u>i  are  the  various  distinct  solutions  of  the  dispersion  relation 
(24)  at  o  =  4/3. 44>46  It  is  strai ghtforward  to  show  that  all  modes  are 
such  that  Re(u)-j)  <  0  except  one  for  which  to  =  0.  Thus,  all  other  modes 
will  decay  in  the  fast  time  scale  and  only  the  neutral  mode  (w  =  0)  will 
persist.  Our  solution  after  a  short  time  will  look  like 

u-  u0+e  C(T )  aCa-Z.  3/.2l  +  0(e2)  (42) 


From  the  given  initial  conditions  it  is  easy  to  determine  the  initial 
condition  for  the  neutral  mode.  Below  the  SCI,  if  the  initial  conditions 

are  such  that  C(0)  >  -  /c/b,  the  solution  will  evolve  to  the  stable  lower 

branch  in  Fig.  (9)  (Region  I),  while  if  C(0)  <  -  /c/b,  C  +  «°  in  finite 
time.  (Fig.  9,  Region  II).  Blow  up  in  finite  time  also  occurs  above  the 
SCL  for  any  C(0).  This  does  not  mean  that  the  actual  solution  blows  up, 
just  that  it  evolves  to  a  final  state  far  away  from  the  two  steady  branches 
shown  in  Fig.  (9),  and  thus,  is  not  accessible  by  perturbation  theory. 

As  can  be  seen  in  (38),  the  blow  up  above  the  SCL  is  described  by  a 
tangent  function,  therefore,  the  growth  rate  we  find  for  this  case  must  be 
appropriately  interpreted.  Note  that  the  linearized  dispersion  relation 
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seems  to  suggest  an  imaginary  exponential  growth  rate  above  the  space 
charge  limit.3*6  In  view  of  our  results,  we  see  that  this  is  actually 
misleading.  Moreover,  we  find  that  even  below  the  SCL  the  stable  steady 
branch  can  be  destabilized  by  sufficiently  large  perturbations. 

Our  results  agree  with  the  linear  theory,  provided  we  consider  the 
limit  where  the  latter  becomes  applicable.  Thus,  we  must  compare  the 
linear  theory  with  (39)  as  t  -*■  +  ®  (near  the  stable  branch)  and  (40)  as 
t  *  -  ®  (near  the  unstable  branch).  To  demonstrate  this  we  set  t  =  et  in 
(39)  and  consider  the  limit  t  +  +  ».  Then 


We  see  that  the  decay  rate  of  the  perturbations  as  t  -*■  +  ®  is  equal  to 

e  „  2.83c  (44) 

a 

which  is  the  same  as  that  found  by  the  linearized  analysis.6  This  quantity 
is  important  as  it  also  determines  an  initial  "growth  rate"  for  the  jump 
instability  above  the  SCL  described  by  (38).  If  we  substitute  the  original 
dimensional  time  variable  into  our  expressions  and  write  the  deviation  of  a 
from  its  value  at  the  SCL  as 
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c 


(a  -  aSCL) 


1/2 


=  (a 


(45) 


we  find  that  the  “growth  rate"  is  given  by 


D 


(a 


4  >1/2  ScE  vo 
V  T  r 


or 


(46) 


(47) 


Of  course,  for  the  expansions  in  (41)  to  be  valid,  we  must  have  C(t)  <  1/e. 
However,  while  C  is  not  too  large  (38)  gives  a  reliable  estimate  for  the 
growth  rate  of  the  instability. 


39 


IV.  MICROWAVE  GENERATION 


Of  the  several  millimeter  sources  that  are  in  various  stages  of 
development  (see  Figure  11),  the  virtual  cathode  oscillator  (Vircator) 
has  a  combination  of  characteristics  which  recommend  it  for  high  fre¬ 
quency  use.  First,  the  frequency  of  the  vircator  is  tunable  by  changing 
the  magnitude  of  an  imposed  axial  magnetic  field,  eliminating  any 
requirement  to  change  the  physical  structure  of  the  device.  A  single 
vircator  will  be  tunable  over  an  order  of  magnitude  in  frequency  (e.g., 

10  GHz  -  100  GHz).  Second,  the  bandwidth  of  the  generator  can  be  narrow 
or  broad  based  on  magnetic  field  shaping.  Third,  because  the  vircator 
functions  above  the  space-charge  limiting  current  for  the  electron  beam, 
given  efficient  operation,  it  should  be  capable  of  much  higher  power  than 
other  microwave  sources.  Finally,  the  lack  of  passive  resonating  struc¬ 
tures  to  produce  the  transmitted  wave  reduces  the  problem  of  field  emis¬ 
sion.  This  also  increases  the  maximum  possible  generator  power. 

Experimentally,  the  virtual  cathode  has  already  proven  itself  to  be 
a  copious  microwave  source.22-24*26*47  Table  4  lists  experiments 
which  have  been  carried  out  to  date.  With  the  exception  of  the  Didenko 
experiment  at  Tomsk,26  the  frequency  spectra  have  all  had  a  broad  band¬ 
width  and  relatively  low  efficiency.  Nevertheless,  even  at  low  effi¬ 
ciency  the  experiment  at  Harry  Diamond  Laboratories47  using  a  foilless 
diode  (I  MV,  30  kA)  produced  3  GW  of  power  in  the  Ku  band.  It  is  one  of 
the  most  powerful  centimeter  wavelength  microwave  sources  available.  The 
experiments  at  the  Naval  Research  Laboratory22  (Mahaffey  et  al.)  and 
Tomsk  will  be  reviewed  in  detail,  because  these  experiments  had  the  most 
useful  diagnostics.  Based  on  analytical  work  and  computer  simulations 
these  experimental  results  can  now  be  explained.  It  will  be  shown  later 
that  the  foilless  diode  in  a  shaped  axial  magnetic  field  represents  the 
best  configuration  for  a  high  frequency  device. 
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Figure  11.  Anticipated  optimal  operating  ranges  for  various  sources  of  high  power 
electromagnetic  radiation.  Several  experimental  points  are  included. 


HIGH  POWER  MICROWAVE  GENERATION  HAS  BEEN  WITNESSED  IN  VARIOUS 
ELECTRON  BEAM  CONFIGURATIONS  WHEN  VIRTUAL  CATHODES  ARE  FORMED. 


DETECTOR  BANDWIDTH  IN  PARANTHESES 


Although  only  scaling  relations  are  presently  available,  qualitative 
dependencies  of  virtual  cathode  parameters  on  beam  kinetic  energy  and 
injected  current  are  known. 


1.  Potential  amplitude,  position,  and  oscillation  frequency  all 
have  the  sane  functional  dependence  on  injected  beam  cur¬ 
rent.  The  positional  dependence  of  the  virtual  cathode  is 
given  in  cigure  12  for  a  50  kV  diode. 


2.  These  parameters  asymptotically  approach  a  limiting  value  for 
current  above  the  space-charge  limit. 

3.  The  fundamental  oscillation  frequency  is  approximately  the 
relativistic  beam  plasma  frequency  given  by 


(48) 


where  n^  is  the  electron  bean  number  density  at  injection.  In  particular 
the  oscillation  frequency  from  numerical  simulations  varies  such  that 
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The  end  result  is  creation  of  an  electrostatic  potential  minimum 
which  oscillates  in  time  and  space.  This  fluctuating  potential  barrier 
acts  as  a  gate  to  reflect  some  electrons  and  transmit  others.  The  motion 
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Variation  of  virtual  cathode  distance,  £,  from  the  one¬ 
dimensional  diode  center  versus  the  difference  between  the 
injected  current,  60,  and  the  current  needed  for  bifurcation 
fig.  ^  units  of  the  Chlld-Langmuir  value.  6b  is  assumed  to 
have  a  value  of  4  from  the  classical  analysis. 


of  the  gate  bunches  charge.  However,  the  charge  bunch  and  virtual  cath¬ 
ode  are  separated  spatially  as  shown  in  Figure  13. 

By  analogy  this  configuration  represents  an  LC  oscillator.  The 
virtual  cathode  acts  as  a  capacitor  to  store  the  beam  kinetic  energy. 
During  that  portion  of  the  limit  cycle  in  which  the  potential  is  greater 
than  the  injected  bean  energy,  charge  is  constrained  to  remain  near  the 
anode.  This  starves  the  virtual  cathode  so  that  its  amplitude  decreases 
below  (yc  _i)  mc2/e.  Once  this  occurs  the  charge  bunch  is  trans¬ 
mitted.  The  electron  motion  represents  a  large  time  varying  current 
through  an  inductor.  The  presence  of  charge  away  from  the  anode  reestab¬ 
lishes  the  virtual  cathode,  and  the  cycle  repeats.  The  effect  on  beam 
current  can  be  examined  using  simulations.  A  net  current  diagnostic  is 
given  in  Figure  14  where  the  probe  is  positioned  between  the  anode  and 
virtual  cathode.  The  injected  current  in  units  of  mc3/e,  v0,  is  three 
times  the  limiting  current,  ,  and  y0  is  5.  Note  that  the  virtual 
cathode  can  actually  reverse  the  direction  of  current.  The  average  cur¬ 
rent  value  is  vA. 

The  oscillating  current  generates  microwaves.25  The  wave  frequency 
is  the  oscillation  frequency  of  the  virtual  cathode.  The  wave  propagates 
down  the  drift  tube  in  a  TM  waveguide  mode,  which  determines  the  wave¬ 
length  and  phase  velocity  of  the  wave  in  the  guide.  The  field  configura¬ 
tion  is  evident  in  simulations  where  there  is  no  axial  magnetic  field. 

If  a  cold  beam  is  injected  and  azimuthal  symmetry  is  assumed  by  the  code, 
the  only  nonzero  fields  are  Ez,  Er,  and  B0  in  cylindrical  geom¬ 
etry.  These  three  fields  define  a  TM  wave  traveling  in  the  z  direction. 

A  time  history  and  frequency  spectrum  of  the  Be  field  downstream  of  the 
virtual  cathode  is  given  in  Figure  15.  The  oscillation  frequency  of  the 
virtual  cathode  and  its  first  three  harmonics  are  clearly  visible.  Also 
included  is  a  history  of  the  Ez  field  on  axis  later  in  time.  The 
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Figure  1 


CCUBE  diagnostic  of  charge  density,  05,  and  potential, 
as  a  function  of  a  longitudinal  distance,  z,  for  an  elec¬ 
tron  beam  in  a  two-dimensional  drift  spcce  with  an 
imposed  large  axial  magnetic  field.  05  is  normalized  to 
its  value  at  injection.  The  value  of  (v0-l)  is  4.  Note 
that  the  minimum  charge  density  and  potential  are 
separated  in  z. 
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Figure  14.  Net  current  amplitude  and  spectrum  for  a  probe  placed  between  t 
anode  and  virtual  cathode,  y  =  5  and  vq  =  3  v?. 


Figure  15.  Be  and  Ez  field  probe  histories  and  frequency  spectrum  for  a  simu 
lation  with  y0  «  5,  v0  *  3.7  and  Bz  *  0.  Field  units  are  MeV/cm 
for  a  beam  density  of  1012. 


frequency  spectrum,  which  is  not  included,  is  almost  identical  to  that  of 
Bq.  No  attenuation  of  the  field  strengths  is  noted  over  the  length  of 
the  drift  tube  which  is  equivalent  to  one  meter  tor  a  beam  density  of 
1012  cm-3.  The  waveguide  mode  is  TMq2  with  an  RF  conversion  efficiency 
of  greater  than  20%.  Because  there  is  no  axial  magnetic  field  in  this 
case  few  electrons  are  transmitted  to  the  end  of  the  drift  tube.  The 
loss  of  transmitted  electrons  does  not  reduce  the  RF  efficiency. 

As  noted  earlier,  several  experiments  have  already  generated  high 
power  microwaves  from  virtual  cathodes.  The  most  informative  experiments 
took  place  at  the  Naval  Research  Laboratory  and  at  the  Institute  of  Nuc¬ 
lear  Physics  at  Tomsk.  Both  experiments  used  a  reflex  triode  configura¬ 
tion.  The  electron  beam  machines  were  almost  identical  (350  kV,  50  kA 
versus  450  kV,  65  kA).  Some  salient  features  were 

1.  The  average  microwave  frequency  increased  as  the  square  root 
of  triode  voltage,  V  (see  Figure  16). 

2,.  Enlarging  the  anode  cathode  gap  spacing,  d,  decreased  wave 
frequency. 

3.  The  microwave  bandwidth  is  related  to  the  change  in  voltage 

over  the  beam  pulse,  aV/V.  In  the  Russian  experiment  where 
A V/V  -  5%  -  10%  the  power  FWHM  occurred  at  ±10%  of  the 

peak  frequency. 

4.  For  a  constant  voltage  the  frequency  decreased  as  cathode 
shape  was  changed  from  solid  to  thick  annular  and  finally 
thin  annular  plates  of  comparable  area. 
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5.  The  imposition  of  an  axial  magnetic  field  greatly  decreased 
the  RF  signal  detected. 

Each  result  can  be  explained  in  terms  of  virtual  cathode  character¬ 
istics.  For  a  nonrelati vistic  beam  in  a  semi -i nfinite  one-dimensional 
drift  space  the  vircator  frequency,  a)QSC,  is  directly  proportional  to 
the  beam  plasma  frequency  given  by  Equation  48.  This  is  a  good  approxi¬ 
mation  for  both  experiments  being  considered  here,  because  the  beams  are 
of  low  energy  and  the  drift  tube  length  is  much  greater  than  its  radius. 
Therefore,  u>osc  a  n^1/2.  The  current  density  extracted  from  the 
cathode  is  defined  by  the  nonrelativistic  Chi  1 d-Langmui r  relation,  which 
is 
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4ird^ 
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An  expression  for  velocity  is  obtained  from  conservation  of  energy.  It 
i  s 


V  ■  (*f  »1/2 

Combining  Equations  3  and  4  and  solving  for  n^,  one  obtains 
1  V 

nb  = 

d 


(51) 


(52) 


or 
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(53) 
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Thus,  cj0SC  a  (V/d2)1/2,  which  is  precisely  the  experimental  result. 

In  addition,  this  dependence  of  vircator  frequency  on  voltage  explains 
the  importance  of  a  flat  top  voltage  pulse  in  producing  a  narrow  band¬ 
width.  The  closure  of  the  A-K  gap,  which  is  due  to  plasma  motion  inher¬ 
ent  in  any  foil  system,  also  makes  the  spectrum  broader.  The  Tomsk 
experiment  used  a  mesh  anode  which  would  minimize  the  latter  effect.  An 
expression  for  these  observations,  which  is  obtained  by  differentiating 
Equation  53,  is 


Ago 
J0SC 


OSC  aV  ,  Ad 

77^  =W  +  T 


(54) 


Cathode  shape  plays  a  role  in  defining  the  space-charge  limiting 
current  via  a  two-dimensional  geometrical  factor.  The  value  of  for 
a  solid  bean  is  less  than  for  a  thin  annular  bean  of  the  same  energy, 
current,  and  area.  In  the  solid  beam  case  the  position  of  the  virtual 
cathode  is  closer  to  the  anode.  This  leads  to  a  higher  oscillation  fre¬ 
quency,  because  the  electron  transit  time  between  the  anode  and  virtual 
cathode  is  less. 

The  final  characteristic  of  the  experiments  is  the  apparent  reduct¬ 
ion  of  microwave  generation  in  the  presence  of  an  axial  magnetic  field. 
However,  two  dimensional  simulations,  which  include  the  diode  physics, 
show  that  microwaves  are  still  being  produced.48  The  magnetic  field 
forces  the  reflected  electrons  to  reenter  the  A-K  gap  region  along  field 
lines.  In  the  case  of  a  foil  diode  or  reflex  triode  the  electrons  suffer 
energy  losses  and  are  scattered  by  the  foil.  After  several  foil  transits 
the  beam  electron  distribution  is  nearly  Maxwellian.  This  is  equivalent 
to  having  electrons  of  widely  varying  voltage  interacting  with  the  virt¬ 
ual  cathode.  Thus,  an  extremely  broad  RF  spectrum  is  produced  which 
greatly  reduces  the  detectable  signal  in  any  particular  frequency  band. 
The  same  result  also  occurs  in  foilless  diodes.  The  broad  electron 
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velocity  distribution  results  from  a  two  stream  instability  between  part¬ 
icles  with  positive  and  negative  axial  velocities  in  the  region  between 
the  real  and  virtual  cathodes.  The  instability  heats  the  two  beams  until 
there  is  a  homogeneous  electron  distribution  in  velocity  space. 

Based  on  our  current  theoretical  knowledge  of  the  virtual  cathode, 
numerical  simulations,  and  experimental  results,  the  following  points 
must  be  considered  in  the  design  of  a  coherent,  high  frequency  vircator. 

1.  Both  the  oscillation  frequency  and  net  current  asymptotically 
approach  a  value  as  injected  current  is  increased  above  the 
space-charge  limit.  Thus,  using  a  large  value  of  v0/v£ 

does  not  substantially  increase  frequency  or  RF  efficiency 
where  efficiency  is  given  by 

n  <  - ^ - o  (55) 

(y0-1)  mc^ 

A  foilless  diode  in  a  strong  axial  magnetic  field  produces  a 
very  thin  annular  beam.  Since  vA  for  an  annular  beam  is 
larger  than  for  a  solid  beam  of  the  same  area,  the  value  of 
v0/v£  will  be  smaller  for  the  same  beam  current. 

2.  Generation  of  high  microwave  frequencies  requires  large  beam 
densities.  Extremely  high  beam  densities  (n  >  1014  cnr3) 
have  been  obtained  from  a  foilless  diode.48  The  beam  plasma 
frequency  scales  linearly  with  the  electron  cyclotron  fre¬ 
quency  due  to  the  magnetic  field.49 
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3.  The  following  characteristics  must  be  met  to  assure  narrow 
bandwidth  microwave  generation  from  the  vircator. 


a)  No  reflexing  of  electrons  in  the  region  between  the  real 
and  virtual  cathodes  must  occur.  Any  axial  magnetic 
field  must  be  shaped  to  divert  the  electrons,  or  flux 
excluders  must  be  employed  to  confine  the  magnetic  field 
to  the  diode  region.  The  latter  arrangement  will  allow 
the  radial  space-charge  electric  field  to  perform  the 
role  of  expelling  electrons  to  the  waveguide  wall. 

b)  The  electron  beam  must  be  cold.  Experimental 50  and  theo¬ 
retical49  results  indicate  that  foilless  diodes  create 
low  emittance  beams.  Laminar  flow,  where  the  electron 
Larmor  orbit  is  smaller  than  the  beam  thickness,  is 
obtained  when50 
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where  a  is  the  orbit  radius,  6  is  the  radial  spacing 
between  the  cathode  and  drift  tube  wall  (which  acts  as 
the  anode)  and  uic  is  the  electron  cyclotron  frequency 
given  by  eBz/mc.  Low  beam  scatter  is  also  assured, 
because  of  the  lack  of  a  foil. 


c)  The  diode  voltage  and  injected  current  must  be  constant. 
Flat-top  voltage  pulses  can  be  attained  in  a  variety  of 
ways  in  several  diode  configurations.  However,  absence 
of  diode  closure  in  some  foil  less  diode  experiments  makes 
a  long  pulse  device  possible.  Also,  at  high  voltages  the 
foilless  diode  operates  as  a  purely  resistive  load. 


therefore 

osc 


constant. 
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It  is  evident  from  this  discussion  that  the  foil  less  diode  in  a 
strong  axial  magnetic  field  represents  the  optimal  configuration  for  a 
high  frequency  vircator.  It  optimizes  microwave  power  and  efficiency 
while  generating  high  frequency,  coherent  radiation. 
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V.  Experimental  Apparatus 


1-  Introduction 


This  section  decribes  the  VIRCATOR  experimental  apparatus  in  detail. 
The  general  design  philosophy  has  been  to  include  as  much  flexibility  as 
possible  in  the  mechanical  and  electrical  systems.  This  approach  allows 
for  the  inevitable  changes  in  experimental  geometry  and  diagnostics  to  be 
made  easily,  and  will  also  permit  a  straight-forward  maximization  of  the 
output  microwave  power  over  the  various  degrees  of  freedom  available. 

Mechanical  aspects  of  the  experiment  are  described  in  Section  2. 
Included  are  the  vacuum  vessel  design  and  descriptions  of  the  anode  posi¬ 
tioning  assembly  and  the  cathode.  Electrical  systems  are  summarized  in 
Section  3.  Electronic  schematics  for  the  pulsed  50  kilogauss  magnet  power 
supply  and  for  the  lumped  element  Blumlein  electron  beam  driver  are  pre¬ 
sented.  Diagnostics  of  the  pulsed  power  systems  and  of  the  electron  beam 
are  also  described.  Measurements  of  the  microwave  emission  power  spectrum 
will  be  made  with  the  microwave  spectrometer  described  in  Appendix  A. 


2.  Mechanical  Design 


The  successful  operation  of  the  VIRCATOR  device  requires  that  the 
electrostatic  potential  surfaces  have  a  particular  structure,  and  that 
these  equipotential  surfaces  be  appropriately  oriented  with  respect  to  the 
applied  magnetic  field.  These  criteria  will  be  satisfied  if  the  cathode, 
anode,  annular  limiter,  and  solenoidal  magnetic  field  are  aligned  concen¬ 
trically,  and  if  all  are  cyl i ndri cal ly  symmetric.  However,  the  required 
accuracy  of  this  concentric  relationship  has  yet  to  be  thoroughly  investi¬ 
gated,  since  it  requires  making  several  runs  of  a  fully  three-dimensional 
simulation  code.  Rather  than  attempt  to  fabricate  the  entire  experiment  to 
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exacting  tolerances,  it  was  therefore  decided  to  build  a  mechanical  system 
capable  of  changing  the  geometrical  relationships  between  the  various  com¬ 
ponents  in  small  incremental  steps.  This  approach  allows  performance  of 
the  experiment  to  be  optimized  in  a  systematic  fashion,  and  will  also  pro¬ 
vide  a  measurement  of  the  system  sensitivity  to  geometrical  relationships 
between  the  components. 

Vacuum  System 


The  components  associated  with  electron  beam  production  are  enclosed 
in  a  standard  18  inch  diameter  bell  jar  vacuum  system,  and  pumped  by  a  NRC 
4  inch  diffusion  pump  to  a  base  pressure  less  than  10~6  torr.  In  order  to 
allow  for  high  voltage  and  mechanical  vacuum  feedthroughs,  a  custom  stain¬ 
less  steel  collar,  designed  by  MRC  and  fabricated  by  Huntington  Mechanical 
Laboratories,  is  positioned  between  the  baseplate  and  glass  bell  jar. 

Design  drawings  for  this  unit  are  shown  in  Figure  17.  Conflat  type  vacuum 
flanges  are  used  on  all  ports  because  of  the  wide  variety  of  electrical  and 
mechanical  feedthroughs  manufactured  for  this  style  flange.  Use  of  off- 
the-shelf  type  hardware  yields  cost  savings  and  increased  system  flexibil¬ 
ity  r'cl:ti*c  lC  custom  machined  vacuum  feedthroughs. 

The  experiment  is  arranged  on  the  feedthrough  collar  of  Figure  17  as 
follows.  A  135  kV  high  voltage  feedthrough  (Ceramaseal  Inc.  #90801854-1) 
is  mounted  on  vacuum  port  1,  positioning  the  carbon  cathode  at  approxi¬ 
mately  5  cm  from  the  center  of  the  chamber.  The  pulsed  magnet  is  mounted 
to  the  stainless  steel  brackets  shown  so  that  it  is  concentric  with  the 
cylindrical  cathode.  The  brackets  are  slotted  to  allow  for  magnet  posi¬ 
tioning  along  its  axis;  vertical  adjustments  are  made  by  screws  in  the 
magnet  base.  The  beam  target  assembly,  consisting  of  the  stainless  steel 
anode,  annular  limiter,  and  waveguide  entrance  port,  is  located  concentric¬ 
ally  with  the  cathode  and  magnet  by  means  of  a  positioning  assembly 
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FEEDTHROUGH  COLLAR 


Figure  17.  Fabrication  drawing  of  the  VIRCATUR  vacuum  vessel.  All  f 
are  Conflat  type  to  accept  standard  feedthrough  hardware. 


FEEDTHROUGH  COLLAR 


gure  17.  Fabrication  drawing  of  the  VIRCATUK  vacuum  vessel.  All 
are  Conflat  type  to  accept  standard  feedthrough  hardware 


FEEDTHROUGH  COLLAR 


Concluded 


mounted  in  vacuum  port  3.  Microwave  emission  into  the  microwave  spectrom¬ 
eter  is  through  vacuum  port  5.  The  remaining  vacuum  ports  accomodate  the 
electrical  and  mechanical  feedthroughs  related  to  beam  and  microwave  diag- 
nosti  cs. 

Beam  Target  Assembly 


The  ideal  experimental  conf i guration  would  allow  manipulation  of  the 
various  electrical  components  from  outside  the  vacuum  vessel,  eliminating 
the  need  to  break  vacuum  when  an  adjustment  is  required.  Because  mechan¬ 
ical  motion  vacuum  feedthroughs  are  costly,  the  design  allows  only  the  beam 
target  assembly,  consisting  of  the  anode,  annular  limiter,  and  the  wave¬ 
guide  entrance  port,  to  be  manipulated  from  outside  the  vacuum  vessel.  A 
commercial  micrometer  positioning  assembly  capable  of  100  micron  spatial 
resolution  (Huntington  Mechanical  Laboratories  #PM-600-XYZTR)  has  been 
chosen  for  this  application.  This  unit  allows  full  5  degree-of-f reedom 
motion  of  the  target  assembly  relative  to  the  catnode. 

The  relative  geometry  of  the  cathode,  beam  target  assembly,  and  magnet 
is  shown  in  Figure  18.  The  cathode  consists  of  a  hollow  carbon  rod  of  out¬ 
side  diameter  0.70  cm  with  a  0.08  cm  diameter  wall.  Two  chacteri sti cs  make 
carbon  a  strong  choice  for  the  cathode.  First,  its  high  resistivity  allows 
magnetic  field  penetration  on  a  sub-microsecond  time  scale,  ensuring  that 
the  magnetic  field  lines  will  be  nearly  normal  to  the  cathode  surface  at 
the  time  that  the  beam  fires.  Second,  it  has  a  low  macroscopic  electric 
field  threshold  for  field  emission,  of  order  10  kV/cm51.  However,  for  the 
design  operating  voltage  of  50  kV,  the  electric  field  at  the  cathode  sur¬ 
face  is  625  kV/cm,  ensuring  that  cathode  field  emission  can  be  obtained  at 
substanti al ly  lower  voltages  and  promoting  uniform  field  emijSion.  The 
purpose  of  the  thin  cathode  wall  is  to  promote  formation  of  an  annular 
electron  beam  by  eliminating  emission  from  the  cathode  center. 
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BEAM  TARGET  ASSEMBLY 


Sketch  of  the  beam  target  assembly.  The  anode,  annular  limiter, 
and  microwave  entrance  port  are  operated  at  ground  potential, 
and  the  c)thode  is  opera' ed  at  negative  high  voltage. 


In  order  to  maximize  the  space-charge  limited  current  emission  from 
the  cathode,  the  cylindrical  anode  of  inside  diameter  0.43  cm  surrounds  the 
cathode  with  a  narrow  0.08  cm  gap.  This  geometry  yields  a  diode  impedance 
for  a  50  kV  beam  of  66  ohms,  with  a  space-charge  limited  current  of  760 
amps.  The  annular  limiter  serves  a  dual  role.  It  acts  as  ground  plane  for 
electron  beam  production  in  the  source  region  and  also  prevents  electrons 
reflected  from  the  virtual  cathode  from  re-entering  the  source  region. 

Beam  equilibrium  considerations  dictate  an  axial  spacing  between  the  annu¬ 
lar  limiter  and  cathode  of  several  anode  radii,  yielding  a  design  gap  of 
2.0  cm.  The  radial  gap  of  the  annular  limiter  is  chosen  to  match  the 
expected  electron  beam  thickness,  0.035  cm,  and  is  centered  just  outside 
the  cathode  radius  at  0.37  cm.  The  diameter  of  the  cylindrical  waveguide 
entrance  port  in  the  virtual  cathode  region  is  chosen  as  1.06  cm,  causing 
the  space-charge  limited  current  in  this  region  to  drop  to  45%  of  its  value 
in  the  source  region.  The  ensuing  virtual  cathode  formation  drives  the 
production  of  50  GHz,  6  mm  wavelength  microwaves  for  a  30  xiiogauss  applied 
field. 

The  components  of  the  beam  target  assembly  of  Figure  18  are  all  nomin¬ 
ally  grounded,  with  the  cathode  driven  at  negative  high  voltage.  In  order 
to  diagnose  the  electron  beam  operation,  however,  the  anode,  annular  limi¬ 
ter,  and  waveguide  entrance  port  are  mounted  in  such  a  way  that,  although 
rigidly  attached  as  a  unit,  the  components  are  electrically  floating.  In 
this  way  the  current  to  each  of  the  electrodes  can  be  monitored,  either  by 
use  of  a  shunt  or  with  a  Rogowski  coil.  Comparing  the  cathode  current  to 
the  anode  current  will  then  yield  an  estimate  of  shank  emission  and  gap 
closure  time;  the  current  to  the  annular  limiter  and  waveguide  entrance 
port  will  give  an  estimate  of  how  well  the  assembly  is  centered  as  well  as 
the  beam  thickness.  These  simple  beam  diagnostics  will  be  supplemented  by 
a  Faraday  cup  current  collector  mounted  on  a  mechanical  probe  drive  mechan¬ 
ism  for  making  detailed  measurements  of  the  electron  beam  profile. 
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Magnet  Design 


The  main  design  criteria  for  the  VIRCATOR  magnet  are  related  to  field 
intensity  and  field  homogeneity.  The  maximum  design  field  is  60  kilogauss, 
required  to  drive  100  GHz  microwaves  in  the  above  experimental  configura¬ 
tion.  Since  diverging  magnetic  field  lines  increase  beam  divergence,  mag¬ 
netic  field  homogeneity  is  require^  in  order  that  the  beam  particle  density 
at  injection  into  the  virtual  cathode  region  be  the  same  as  the  beam  elec¬ 
tron  density  it  the  source.  An  additional  constraint  imposed  on  the  design 
is  that  the  working  magnetic  field  volume  be  several  times  larger  in  radius 
and  length  than  the  nominal  anode  size  described  in  the  previous  section. 
The  magnet  has  been  designed  oversize  to  allow  for  the  contingency  that 
anode-cathode  gap  closure  problems  with  the  nominal  0.08  cm  wide  gap  may 
require  using  a  much  larger  anode-cathode  radial  spacing,  of  order  milli¬ 
meters.  Maintaining  diode  impedance  in  the  50  ohm  range  would  then  require 
increasing  the  anode  diameter  to  centimeters.  The  present  magnet  design 
admits  this  option. 

The  magnet  overall  length  is  15.2  cm  with  a  usable  inside  diameter  of 
3.5  cm.  To  produce  a  60  kilogauss  field  on  axis  at  the  center  of  the  mag¬ 
net  requires  a  linear  current  density  of  49,500  amp-turns/cm.  The  magnet 
is  wound  with  284  turns  of  1  mm  diameter  solid  copper  wire  so  that  the  coil 
current  is  2.6  kA  and  magnet  inductance  is  1  millihenry.  In  order  to 
enhance  field  homogeneity,  the  coil  is  wound  with  additional  turns  at  the 
ends.  The  self  forces  on  the  magnet  tend  to  squeeze  it  axially  and  expand 
it  radially.  Sel f-compressi on  forces  in  the  axial  direction  are  taken  up 
by  the  copper  windings  and  the  epoxy  binder,  and  the  radial  expansion 
forces  are  born  by  an  encapsulating  cylindrical  stainless  steel  tube. 
Neglecting  the  restoring  forces  provided  by  the  copper  windings  and  epoxy 
binder  yields  a  stress  of  15,300  psi  in  the  3.2  mm  thick,  6.4  cm  outside 
Hiameter  tube  for  a  60  kG  field.  Since  the  yield  strength  of  type  304 
stainless  steel  exceeds  40,000  psi,  the  design  provides  an  acceptable 
safety  factor. 
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The  magnetic  field  risetime  is  determined  by  the  magnetic  inductance, 
source  capacitance,  and  series  resistance.  As  described  in  the  next 
section,  these  were  chosen  to  give  a  300  microsecond  risetime.  Since  the 
field  penetration  time  for  the  anode,  cathode,  and  the  stainless  steel 
magnet  casing  are  all  less  than  15  microseconds,  the  magnetic  field  in  the 
electron  beam  source  region  can  be  accurately  calculated  using  computer 
codes  which  neglect  eddy-current  effects. 

3.  Electrical  Design 


The  VIRCATOR  electrical  system  consists  of  the  magnet  capacitive  dis¬ 
charge  circuit,  the  lumped. element  Blumlein  pulse  forming  line  for  driving 
the  electron  beam,  and  associated  timing  and  data  recording  channels.  In 
order  to  minimize  the  emission  of  electromagnetic  radiation  from  the  high 
speed  switching  circuits,  all  high  voltage  components  are  enclosed  in  a 
grounded  Faraday  cage  formed  from  bronze  screen.  Additional  shielding  for 
the  data  channels,  consisting  of  five  single  trace  oscilloscopes  and 
recording  cameras,  is  provided  by  a  fully  screened  instrumentation  rack. 
Because  of  the  limited  number  of  data  channels,  the  generation  of  a  com¬ 
plete  data  set  will  require  multiple  shots.  Due  to  the  expected  high  shot 
frequency,  one  every  5  minutes,  data  acquisiton  will  still  be  relatively 
rapid,  provided  shot  reproducibility  is  acceptable. 

Magnet  Pulsed  Power  Supply 

A  schematic  for  the  circuit  which  drives  the  magnetic  field  coil  is 
shown  in  Figure  19.  A  90  microfarad  energy  storage  capacitor  is  discharged 
into  the  magnetic  field  coil  through  a  series  resistance  which  is  chosen  so 
that  the  circuit  is  critically  damped.  Since  the  current  does  not  ring, 
the  lifetime  of  both  the  capacitors  and  the  magnetic  field  coil  are 
increased.  The  price  paid  is  an  increase  by  a  factor  of  approximately  2.5 
in  the  required  bank  capacitance  relative  to  a  lightly  damped  circuit. 
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MAGNET  PULSE  POWER  SUPPLY 


Schematic  of  the  pulse  power  circuit  which  drives  the  magnetic 


The  same  result  could  be  achieved  by  crowbaring  the  current  at  peak  field 
with  an  ignitron,  and  eliminating  the  series  resistance.  Since  the  cost  of 
the  increased  capacitance  is  small  compared  to  the  cost  of  the  ignitron  and 
its  driving  circuit,  the  former  option  was  chosen. 

The  magnet  energizing  circuit  of  Figure  19  operates  as  follows.  A  20 
kv,  10  mA  power  supply  (Del  Electronics  #R IU-20-10)  charges  the  90  micro¬ 
farad  capacitor  bank  to  its  charge  voltage,  about  10  kv  for  a  40  kilogauss 
field.  A  krytron  gas  tube  switch  is  command  triggered  to  discharge  a  1 
microfarad  capacitor  charged  to  2.5  kv  into  the  ignitor  terminals  of  a 
mercury  ignitron  (General  Electric  #GL-7703).  The  triggered  ignitron  then 
discharges  the  main  capacitor  bank  into  the  inductive  load  through  a  6.7 
ohm  series  resistor.  The  peak  magnetic  field  occurs  at  300  microseconds, 
and  decays  with  a  time  constant  of  300  microseconds.  A  time  delay  gener¬ 
ator  (California  Avionics  #1 13CR )  command  fires  a  spark  gap,  initiating  the 
electron  beam  at  peak  magnetic  field  time. 

Electron  Beam  Pulsed  Power  Supply 


The  ideal  voltage  pulse  shape  for  driving  the  electron  beam  would  have 
a  fast  leading  edge  for  initiating  field  emission  at  the  cathode,  followed 
by  a  somewhat  lower  voltage  flat  top  for  driving  constant  current  in  an 
assumed  constant  impedance  load  provided  by  the  diode.  A  good  apporoxima- 
tion  to  this  ideal  source  is  provided  by  a  lumped  element  artificial 
Blumlein  pulse-forming  line  shown  in  Figure  2G,  which  uses  high  energy 
density  barium  titanate  capacitors  as  the  energy  storage  medium.  This  cir¬ 
cuit  provides  twice  the  charge  voltage  into  an  open  load,  thereby  initiat¬ 
ing  breakdown  at  the  cathode,  followed  by  an  approximately  flat  top  signal 
equal  to  the  charge  voltage  when  driving  a  matched  load.  For  a  N-stage 
pulse-line  the  risetime  x,  pulse  length  T,  and  impedance  Z  are  given  by 
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ELECTRON  BEAM  PULSED  POWER  SUPPLY 


chemacic  of  the  lumped  element  Olumlein  pulse-forming  line, 
nitally  10  stages  will  be  used  to  give  a  pulse  length  of  600 
anoseconds . 


t  =  /nr? 

(57) 

T  =  2N/TT 

(58) 

Z  =  2  STJZ 

(59) 

where  L  and  C  are  the  component  inductances  and  capacitances,  respec¬ 
tively.  In  order  to  build  in  future  capability  for  100  kV  operation,  three 
40  kV  disc  capacitors  will  be  connected  in  series  to  give  900  pF  rated  at 
120  kV.  Matching  the  pulse-line  impedance  to  the  desired  diode  impedance 
of  66  ohms  requires  an  inductance  of  .98  microhenries.  For  a  ten  stage 
pulse  forming  line  the  risetime  and  pulse  length  follow  from  above  as  15 
and  594  nanoseconds,  respectively.  Increasing  the  number  of  stages  to  pro¬ 
vide  full  1  microsecond  pulse  length  capability  is  straight-forward  and 
will  be  done  after  any  gap  closure  problems  occuring  at  the  shorter  pulse 
lengths  have  been  solved.  Increasing  the  number  of  stages  has  an  addi¬ 
tional  benefit  of  smoothing  the  flat  top  portion  of  the  voltage  pu1?®. 

The  trigger  and  charging  circuit  for  the  electron  beam  driver  are 
shown  in  Figure  21.  The  Blumlein  array  is  charged  by  a  90  kV  power  supply 
(Del  Electronics  #R IU-90- 1- 1 )  and  command  fired  at  the  time  of  peak  mag¬ 
netic  field  by  a  signal  from  the  time  delay  generator.  The  time  delay  gen¬ 
erator  triggers  a  55  kV  pulse  generator  (Pacific  Atlantic  Electronics  #PT- 
55-M)  which  in  turn  drives  the  spark  gap  (Pulsar  #670).  The  spark  gap  dis¬ 
charges  the  pulse-forming  line  into  the  load  provided  by  the  electron 
beam.  A  Rogowski  coil  and  voltage  divider  monitor  the  diode  current  and 
voltage,  respectively. 
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VI.  TECHNICAL  REQUIREMENTS  AND  TASKS 

HOW  THE  STATEMENT  OF  WORK  HAS  BEEN  SATISFIED 

In  this  section  we  review  the  tasks  encompassed  by  the  statement  of 
work  and  detail  how  these  technical  requirements  have  been  satisfied.  We 
also  include  a  list  of  publications  and  talks  given  on  this  subject,  pro¬ 
fessional  personnel  associated  with  the  research  effort,  and  consultative 
and  advisory  functions. 

Statement  of  Work  (0001AA) 

1.  Theoretically  determine  an  optimal  configuration  for  a  proof -of- 
principal  vircator  which  operates  in  the  millimeter  wavelength  regime. 

We  have  reviewed  in  detail  how  one  can  accurately  calculate  the  space- 
charge  limiting  current,  which  is  a  necessary  requirement  for  determining 
if  a  virtual  cathode  will  form,  in  Section  II.  For  the  first,  time  both  the 
nonlinear  treatment  of  the  space-charge  limit  instability  and  the  scaling 
of  key  parameters  (oscillation  frequency,  potential  depth,  location  of  the 
virtual  cathode)  has  been  accomplished.  This  is  presented  in  Section  III. 
Work  is  now  in  progress  cn  the  effect  of  electron  temperature  and  multiple 
virtual  cathode  formation.  (See  list  of  proposed  publications.) 

It  is  now  clear,  based  on  simulations  performed  at  Los  Alamos  National 
Laboratory  (LANL)  and  from  our  analytical  efforts  that  beam  temperature  not 
only  effects  virtual  cathode  frequency  bandwidth,  but  more  importantly 
microwave  generation  efficiency.  The  conclusion  of  this  report  that  a 
toilless  diode  and  some  means  of  preventing  electron  reflexing  are  optimal 
for  producing  a  narrowband,  efficient  oscillator  is  certainly  justified. 
This  work  is  presented  in  Section  IV. 
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2.  Accomplish  the  technical  design  of  the  proof-of-pri nciple  vir- 
cator. 

The  technical  design  of  the  vircator  is  given  in  Section  V.  Of  great¬ 
est  importance  is  the  flexibility  which  is  provided  by  this  plan.  For 
instance,  the  field  coil  is  devised  to  allow  a  60  kG  field  with  an  accept¬ 
able  margin  of  mechanical  safety.  Also  the  artificial  Blumlein  willallow 
us  to  obtain  100  kV  pulses,  if  we  used  three  rather  than  two  barium  titan- 
ate  capacitors  per  stage. 

3.  Fabricate  portions  of  the  vircator  incorporating  long  lead 
items. 

We  have  designed  the  custom  bell  jar  vacuum  feedthrough  system.  It  is 
currently  being  manufactured.  All  of  the  electrical  components  (capaci¬ 
tors,  vacuum  feedthroughs,  delay  generators,  switches,  and  power  supplies) 
have  been  specified,  ordered  and  delivered.  Mechanical  components  for 
accurate  positioning  of  the  beam-target  assembly  and  raw  materials  for  the 
assembly  have  been  ordered  and  delivered.  The  construction  of  the  Blumlein 
and  assembly  will  be  completed  shortly. 

4.  Design,  construct  and  calibrate  a  quasi-optical  microwave  spec¬ 
trometer. 

The  design  of  the  spectrometer  is  presented  in  Appendix  A.  It  has 
been  constructed.  Calibration  using  high  frequency  sources  available  at 
the  Naval  Research  Laboratory  (NRL)  will  be  completed  in  the  month  of 
February.  In  order  to  maximize  the  utility  of  the  spectrometer  its  current 
frequency  range  (75-110  GHz)  will  be  extended  down  to  35  GHz.  This  addi¬ 
tion  will  be  accomplished  by  April  1,  1983. 
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5.  Characterize  the  microwave  emission  (power,  frequency  and  band¬ 
width)  of  a  relativistic  electron  beam  in  a  vircator  configuration  using 
the  microwave  spectrometer. 

This  last  task  has  not  been  accomplished.  The  lack  of  availability  of 
time  on  such  a  beam  machine  at  the  Air  Force  Weapons  Laboratory  (AFWL)  and 
the  longer  than  anticipated  time  required  to  complete  task  4,  both  contrib¬ 
uted  to  this  situation.  We  have  postponed  this  requirement  to  CY  1983. 
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Bibliography  of  Talks  and  Publications 
Archival  Scientific  Publications 


"Design  of  a  Millimeter  Wave  Source  Using  a  Virtual  Cathode  Oscillator 
(VIRCATOR)",  in  Conference  Record  -  1982  IEEE  International  Conference  on 
Plasma  Science,  p.  64. 

"Effects  of  Temperature  on  the  Space-Charge  Limit  of  Electron  Beams",  to  be 
submitted  to  Applied  Physics  Letters  (E.  A.  Coutsias  and  D.  J.  Sullivan). 

"Formation  of  Singularities  in  One-Dimensional  Electron  Flows",  to  be  sub¬ 
mitted  to  Physics  of  Fluids  (E.  A.  Coutsias  and  D.  J.  Sullivan). 

"High  Power  Microwave  Generation  from  a  Vi rtual  Cathode  Oscillator",  to  be 
published  in  IEEE  Trans.  Nuc.  Sci.,  June  1983  (D.  J.  Sullivan). 

Classified  Scientific  Reports 


"Theory  and  Design  of  a  High  Power  Virtual  Cathode  Oscillator  (VI-RCATOR)", 
to  be  published  in  Proceedings  of  the  2nd  National  Conference  on  High  Power 
Microwave  Technology  (D.  J.  Sullivan). 

Internal  MRC  Reports 


"Theory  and  Application  of  the  Virtual  Cathode  in  Charged  Particle  Beams," 
AMRC-R-433  (1982). 

Talks  and  Meetings 

1982  IEEE  International  Conference  on  Plasma  Science,  Ottawa,  Canada 

"Design  of  a  Millimeter  Wave  Source  Using  a  Virtual  Cathode  Oscillator 
(VIRCATOR)" 

MGX  Review  Meeting  -  Lawrence  Livermore  National  Laboratory  (22  Sep  82) 
"MGX  Simulations  Using  CCUBE" 

Virtual  Cathode  Microwave  Design  Review  -  Los  Alamos  National  Laboratory 
(31  Jan  83) 

"VIRCATOR  Theory  and  Design  at  Mission  Research  Corporation" 

"The  Mission  Research  Corporation  VIRCATOR" 

Second  National  Conference  on  High  Power  Microwave  Technology  -  Harry 
Diamond  Laboratories  (Mar  83) 

"Theory  and  Design  of  a  High  Power  Virtual  Cathode  Oscillator  (VIRCATOR)" 

1983  Particle  Accelerator  Conference  -  Santa  Fe,  New  Mexico 

"Hiah  Power  Microwave  Generation  from  a  Virtual  Cathode  Oscillator 
(VIRCATOR)" 
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Professional  Personnel  (Alphabetical  order) 


W.  Michael  Bollen 

Evangelos  A.  Coutsias  ('JNM  Department  of  Mathematics  and  Statistics) 

Carl  A.  Ekdahl 

Robert  H.  Jackson 

Bret  E.  Simpkins 

Donald  J.  Sullivan 

Donald  E.  Voss 

John  Walsh  (Dartmouth  College,  Department  of  Physics  and  Astronomy) 
Consultative  and  Advisory  Functions 


Consulting  Agreement  with  Lawrence  Livermore  National  Laboratory  (LLNL)  to 
aid  in  design  of  the  Microwave  Generator  Experiment  (MGX)  which  uses  a 
VIRCATOR  as  the  microwave  source.  Points  of  contact  at  LLNL  Urs.  Hriar 
CaDa^an,  Wayne  Hohfer,  Raymond  Scarpetti  and  George  Vogtlm.  Dates  pre¬ 
sent  at  LLNL  25-27  May  1982  and  20-23  September  1982. 


Numerous  informal  discussions  with  Los  Alamos  National  Laboratories  (LANL) 
personnel  in  Groups  R-l  and  X-8  on  using  the  VIRCATOR  as  a  microwave 
source.  Tn  particular  discussions  with  Drs.  Carl  Ekdahl  and  Thomas  Kwan. 
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SECTION  1 
INTRODUCTION 


1. 1  Background 


In  characterizing  electromagnetic  radiation,  the  intensity  and 
wavelength  of  the  radiation  are  used.  A  device  for  measuring  these  two 
quantities  is  the  spectrometer.  Spectrometers  are  a  commonly  used  device 
in  the  UV,  visible  and  IR  region  of  the  spectrum.  A  variety  of  techniques 
for  performing  the  spectrometry  may  be  used:  Fabry-Perot  interferometry, 
Fourier  transform  spectroscopy,  and  diffraction  grating  spectroscopy.  For 
the  microwave  spectrum,  the  Fourier  transform  technique,  commonly  called  a 
spectrum  analyzer,  has  been  used.  The  other  methods  have  not  had 
widespread  use  with  microwaves. 

Spectrometry  for  non-repetiti ve,  short  pulse  length  microwave 
radiation  can  only  be  easily  performed  using  the  grating  spectrometry 
method  mentioned  above,  and  it  has  been  used  with  some  success2 » 3 ^ . 

The  main  difficulty  encountered  has  been  in  calibration  of  the 
spectrometer.  Calibration  is  required  because  the  response  typically  is 
not  flat;  the  signal  output  varies  with  wavelength  for  a  constant 
intensity  input.  Further,  calibration  sources  are  not  easily  available 
except  at  very  low  power  levels  (milliwatts).  In  our  design  we  have 
attempted  to  make  the  response  flat;  however,  we  are  limited  by  the 
nonuniform  response  of  the  microwave  detectors. 

1.2  The  Diffraction  Grating  Spectrometer 


A  basic  description  of  diffraction  grating  spectrometers  can  be 
found  in  many  optics  texts  (see,  for  exa,  ^le,  reference  5).  The  physics 
of  the  diffraction  process  is  described  by 

mA=d(sina+sin8)  (1) 
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where  m  is  the  order  of  the  diffraction  (0,  1.1,  12,  ...),  X  is  the 
wavelength,  d  is  the  spacing  between  the  grating  grooves,  a  and  g  are  the 
angles  of  incidence  and  diffraction  respectively  (see  Figure  1).  The  plus 
sign  is  used  in  equation  1  if  a  and  g  are  the  same  sides  of  the  surface 
normal;  otherwise,  the  minus  sign  is  used.  Figure  1  depicts  the  commonly 
used  echelette  grating  (ninety  degree  apex  angles).  The  echelette  is 
popular  because  it  is  easy  to  fabricate,  and  it  has  been  shown 
experimentally  to  give  the  best  efficiency  for  a  ruled  grating. 

The  dispersion  of  the  grating  (for  fixed  incident  angle,  a)  is 
expressed  as 

dg/dX  =  m/(d  cos  g)  (2) 

where  terms  are  defined  as  before.  The  resolving  power,  R,  is 

R  =  X/AX  =  mN  (3) 


where  N  is  the  number  of  grooves;  however,  there  is  a  fundamental 
1  imitation  on  R, 

R  1  2W/X  (4) 

where  W  is  the  width  of  the  grating  (W  =  Nd).  The  details  of  our  design 
using  these  equations  may  be  found  in  the  next  section. 
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t i on  of  geometry  for  an  tcheJlette  grating. 


SECTION  2 

DESIGN  OF  A  W-BAND  SPECTROMETER 


2. 1  Design  Philosophy 


A  microwave  spectrometer  has  been  designed  to  measure  frequency 
and  bandwidth  for  radiation  from  the  Mission  Research  Corporation  (MRC) 
vircator  experiment.  An  effort  has  been  made  to  improve  upon  the  basic 
design  currently  in  use2»3>4  by  incorporating  new  knowledge  on 
diffraction  gratings6  and  by  the  use  of  innovative  design  concepts. 
Computational  modeling  suggests  3-4  mm  radiation  from  the  vircator  with 
fine  structure  of  2  GHz.  These  were  used  as  the  basic  design 
constraints.  Additionally,  we  imposed  the  following  constraints  to  make 
the  device  easier  to  use:  1)  limit  the  wavelength  range  to  remove  high 
order  (m  >  1)  diffraction,  2)  demand  uniform  (flat)  device  response,  3) 
require  wide  device  bandwidth,  4)  require  wide  bandwidth  for  any  single 
input  pulse,  5)  allow  for  modification  to  different  wavelength  bands,  6) 
minimize  unwanted  signal  noise,  and  7)  keep  the  device  small  and 
portable.  Clearly,  all  these  constraints  can  not  be  simultaneously  met. 
and  trade-offs  were  required. 


The  final  design  is  shown  in  Figure  2.  This  spectrometer  is 
designed  to  operate  at  W  band  (2.7-4  mm)  with  a  minimum  of  2  GH-' 
resolution  over  the  entire  band.  The  bandwidth  of  the  device  is  limited 
on  the  long  wavelength  end  by  the  W-band  input  waveguide,  and  on  the  short 
wavelength  end  by  the  low-pass  diffraction  grating.  This  effectively 
eliminates  higher-order  (m  >  1)  diffraction.  Care  was  taken  in  the  design 
of  the  diffraction  gratings  to  choose  a  design  with  a  flat  frequency 
response.  The  main  limitation  encountered  in  obtaining  uniform  response 
was  in  the  microwave  detectors.  These  do  not  have  uniform  response,  and 
their  response  varies  with  each  detector.  Wide  device  bandwidth  is 
obtained  by  the  use  nf  a  rotatable  diffraction  grating.  Wide,  single-shot 
bandwidth  is  obtained  by  using  a  large  reflection  mirror  (allowing  for 
more  than  one  diffraction  angle)  and  use  of  a  multichannel  detection 
array.  We  have  chosen  a  six-channel  detection  array  designed  for  10  GHz 
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full-width  band  response  centered  at  4  mm  (75  GHz).  Modification  to 
different  bands  is  possible  by  replacing  the  two  diffraction  gratings. 
Where  possible,  device  size  was  reduced.  For  example,  folding  of  path 
length  was  done.  Path  length  folding  and  use  of  microwave  absorber  also 
removed  any  direct  line  of  sight  between  the  initial  input  signal  and  the 
diffracted  signal  detectors.  Additionally,  a  collimating  horn/lense 
system  was  used  allowing  a  slight  savings  in  size  over  that  achieved  using 
a  mirror  system.  The  details  concerning  the  design  of  individual 
components  is  found  in  the  next  section. 

2.2  Design  of  the  MRC  W-Band  Spectrometer 

2.2.1  Microwave  Input 


The  microwaves  enter  the  spectrometer  through  a  length  of  W-band 
waveguide  (see  Figure  2).  The  waveguide  performs  the  first  filtering  by 
rejecting  microwave  wavelengths  greater  than  4  mm.  The  input  signal  must 
be  attenuated  to  a  level  of  a  few  hundred  milliwatts  to  prevent  damage  to 
the  microwave  detectors.  Some  attenuation  (10-20  dB)  will  occur  in  the 
spectrometer  itself  (due  to  losses);  however,  the  main  attenuation  must  be 
performed  externally.  This  is  performed  geometrically  by  moving  the  input 
of  the  wavequide  away  from  the  output  of  the  microwave  source;  then 
geometrical,  1/RZ,  attenuation  will  occur.  A  10  MW  signal  will  be 
attenuated  to  100  mW  in  the  waveguide  input  for  a  separation,  R,  of  10 
meters . 


Geometrical  attenuation  assumes  a  uniform  distribution  of  the 
microwaves  over  a  spherical  surface  (actually  half  of  the  surface  for  our 
case).  The  accuracy  of  this  assumption  depends  on  two  factors.  First, 
the  spectrometer  must  be  located  in  the  far  field; 

R  >  02/X 


A 


(5) 


where  D  is  the  diameter  of  the  source  aperture,  X  is  the  wavelength  of  the 
radiation,  and  R  is  the  separation  distance.  For  most  cases  of  interest 
to  us,  this  is  usually  easily  satisfied.  For  example,  for  a  1  cm  wide 
output  aperture  at  4  mm,  the  separation  must  be  2.5  cm.  Second,  a  TEn 
mode  (dominant  mode)  in  circular  waveguide  has  been  assumed.  Other  modes 
can  cause  the  distribution  of  power  to  be  peaked  at  certain  angles.  For 
this  reason,  a  mapping  of  the  angular  power  distribution  should  be 
performed.  If  it  is  not  uniform,  then  the  spectrometer  measurements 
should  be  made  at  several  angles.  It  is  conceivable  that  the  frequency  of 
the  microwaves  from  the  source  is  linked  to  the  mode.  Thus,  a  measurement 
of  the  entire  spectrum  requires  a  sampling  of  all  modes. 

Inside  the  spectrometer  the  input  signal  is  expanded  using  a 
standard  gain,  rectangular  horn.  This  horn  has  a  beam  width  (FWHM)  of 
20°.  A  collimating  lense  is  placed  in  the  far  field  of  the  horn  and  at  a 
distance  such  that  half  of  the  input  signal  is  intercepted.  The  lense  is 
a  standard  hyperbolic  design.7  The  front  surface  is  hyperbolic  in  shape 
and  is  described  mathematically  by 

r  .  <n  -  n  1 _  (6) 

n  cos  0-1 

where  n  is  the  index  of  refraction  and  R,  L,  and  0  are  shown  in  Figure  3. 
The  lense  is  made  from  Rexolite  (n  =  1.6).  Rexolite  has  a  large  index  of 
refraction  and  low  attenuation  for  millimeter  microwaves.8  The  horn/lense 
collimation  technique  is  compact  and  aids  in  reducing  the  overall  size  of 
the  spectrometer. 
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2.2.2  Diffraction  Grating 

Although  the  basic  theory  of  diffraction  gratings  has  already 
been  discussed  in  Section  1.2,  a  few  more  definitions  are  required.  When 
discussing  gratings,  the  typical  mount  used  is  the  Littrow  mount.  For 
this  case,  the  incident  and  diffracted  rays  occur  on  the  same  line  ( e  = 
-a).  Under  these  conditions  equation  (1)  simplifies 

mx  =  2d  sin  a  .  (7) 

Further,  if  the  angle  of  the  grooves,  9  (see  Figure  1),  is  the  same  as  a, 
maximum  diffraction  will  occur  since  the  diffracted  and  reflected  angles 
are  the  same.  This  condition  is  referred  to  as  being  "in  blaze"  and  0  is 
the  blaze  angle.  The  Littrow  mount  is  obviously  not  realizable  in  actual 
practice;  however,  for  TM  radiation  (H  parallel  to  the  grooves),  a 
departure  from  Littrow  still  yields  a  high  and  nearly  uniform  efficiency, 
Drovi hed  that  0.85  _<  X/d  _<  1.5.  This  is  true  up  to  45°  off  Littrow. 
Therefore,  the  gratings  were  designed  using  TM  radiation  and  assuming  a 
Littrow  mount.  An  Echelette  grating  design  was  chosen  and  a  blaze  angle 
of  30°  was  picked.  For  this  blaze  angle,  nearly  uniform  transmission 
occurs  for  .7  _<  x/d  _<  i. 7  (see  Figure  4)6.  From  Figure  4  it  is  clear  that 
the  grating  operates  best  for  1  _<  X/d  £  1.5.  This  fixes  d  for  the 
grating.  The  grating  characteristics  are  shown  in  Table  1.  In  order  to 
use  a  rotating  diffraction  grating,  the  angle  a  -  B  was  fixed.  The 
detectors  and  collecting  mirror  can  then  be  fixed  with  respect  to  the 
incoming  microwaves,  and  the  frequency  ic  scanned  by  rotating  the 
grating.  In  our  design  a  -  3  is  set  equal  to  25°.  This  angle  is  kept  as 
small  as  possible  to  stay  near  Littrow,  but  must  be  made  large  enough  so 
the  reflecting  mirror  does  not  occlude  the  diffraction  grating. 

In  addition  to  the  diffraction  grating  which  disperses  the 
microwaves  for  wavelength  measurements,  another  grating  is  used  in  the 
zeroth  order  (m  =  0)  to  limit  the  wavelength  of  the  radiation  to  the 
dispersing  grating.  It  can  be  shown6  that  a  grating  acts  like  a  mirror  in 
m  =  0  for  sufficiently  large  wavelength  radiation.  The  criteria  is 

d  <  Xmin/(1  +  sin  a)  (8) 
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Efficiency  of  Littrow  mount  Echellete  gratings  with 
blaze  angle  9.  The  solid  line  is  TM  pol ari zati on 
(H  parallel  to  grating  grooves)  and  the  dashed  line 
is  TE  polarization  (E  parallel  tc  grating  grooves. 
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TABLE  1 


F 

X/d 

a 

3 

no 

1.01 

43.65 

18.50 

100 

1.11 

47.2 

22.2 

90 

1.24 

52.0 

27.0 

80 

1.39 

58.0 

33.0 

75 

1.48 

61.85 

36.85 

F  is  frequency  a  is  angle  of  incidence 

X  is  wavelength  3  is  angle  of  diffraction 

d  is  groove  spacing 

refer  to  Figure  1  for  definition  of  geometry 

For  the  grating  d  =  0.27  cm 

6  =  30°  (Note:  The  grating  is  also  in  blaze  for 
90  -  6  =  60°) 

a  and  6  are  computed  based  on  a  -  $  ~  2&q  =  25°. 
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where  Xm-jn  is  the  shortest  wavelength  for  which  the  grating  will  act  as 
a  mirror.  We  chose  a  equal  to  45°  and  Amjn  equal  to  2.2  mm  (140  GHz). 
Using  this  grating  as  a  reflecting  mirror  effectively  removes  the 
possibility  of  higher  order  radiation  appearing  on  our  detectors.  (For 
example,  2  mm  radiation  for  m  =  2  will  appear  the  same  as  4  mm  radiation 
for  m  =  1.)  The  mi rror/grating  then  reflects  all  radiation  with  a 
wavelength  greater  than  2.2  mm  and  diffracts  all  shorter  wavelength 
radiation.  We  have  placed  microwave  absorbers  in  the  proper  position  to 
remove  the  diffracted  short  wavelength  radiation. 

2.2.3  Diffracted  Signal  Detection 


As  mentioned  before,  a  multichannel  detection  system  was  decided 
upon  to  give  broadband  response  on  a  single-shot  basis.  It  was  decided  to 
design  for  a  10  GHz  bandwidth  centered  at  75  GHz  (4  mm).  The  detection 
system  consists  of  six  microwave  detectors.  A  spherical  mirror  is  used  to 
concentrate  the  diffracted  radiation.  The  mirror  and  detector  locations 
were  designed  graphically  using  ray  tracing.  The  mirror  is  55  cm  wide  and 
has  a  200  cm  radius  of  curvature.  The  detectors  are  located  on  a  focal 
line  which  is  approximately  100  cm  from  the  mirror.  The  mirror  itself  is 
made  using  thermoform  plastic.  A  wooden  master  is  machined.  The  plastic 
is  then  formed  to  the  master  and  coated  with  a  silver,  high-conductivity 
paint. 


The  actual  signal  detection  is  performed  using  1N53  diode 
crystals  mounted  in  W-band  waveguide.  Although  these  detectors  are  not  as 
flat  in  response  as  others,  their  response  is  not  sufficiently  poor  to 
justify  the  increased  expense  of  more  uniform  response  detectors.  Expense 
becomes  an  important  factor  due  to  the  potentially  short  lifetime  of  the 
detector.  A  large  microwave  signal  permanently  destroys  the  crystal.  It 
has  been  our  experience  in  high-power  microwave  diagnostics  that  diode 
crystals  are  continually  being  burned  out.  This  is  due  to  the  large  power 
fluctuations  that  occur  in  the  experimental  microwave  sources.  Due  to 
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this  non-uniform  response  each  detector  will  have  to  be  individually 
calibrated  for  frequency  response.  The  output  of  the  detectors  will  be 
available  on  50  ohm  BNC  connectors  located  on  the  outside  of  the 
spectrometer.  If  long  cables  will  be  used,  it  is  recommended  that  line 
driver  amplifiers  also  be  used. 
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